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Abstract 


This paper presents a method for estimating the spectra of water wave disturbances on 
five of the six axes of a stationary, slender body underwater vehicle in an inertia dominated wave 
force regime, both in head seas and in beam seas. Inertia dominated wave forces are typical of 
those encountered by a 21 inch diameter, torpedo shaped underwater vehicle operating in coastal 
waters and sea state 2. Strip theory is used to develop transfer function phase and magnitude 
between surface water waves and the slender body pitch, heave, and surge forces and moment for 
the vehicle in head seas, and for pitch, heave, yaw, and sway forces and moments in beam seas. 
Experiments are conducted which verify this method of transfer function calculation, and 
demonstrate the effects of vehicle forward motion in the head seas case. Using known sea spectra 
and linear time invariant systems theory allows for estimation of the water wave disturbance 
spectra for these forces and moments. 

Application of sliding control techniques are then developed for the underwater vehicle 
longitudinal plane. equations of motion. Computer simulations are used to demonstrate the 
dependence of underwater vehicle depth control upon the pitch control, and adaptive pitch control 
is shown to provide good performance in the presence of substantial parametric uncertainty. 
Pitch disturbance rejection properties of variations of the sliding controller are investigated. Both 
single frequency and stochastic disturbances are used, and the stochastic disturbance is developed 
using the results of the earlier investigation. 
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Chapter1 INTRODUCTION 


1.1 Motivation 


Autonomous Underwater Vehicles (AUVs) have become increasingly useful tools in the 
exploration of the ocean depths, where the effects of surface waves are far removed from the 
operating region of the vehicle. As the range of missions for AUVs expands, so does the need to 
understand the disturbances which the vehicle will encounter in its enlarged theater of operation. 
While deep underwater, ocean currents may be the source of the predominant disturbance to the 
untethered AUV, the effect of gravity water waves becomes important when operating an AUV 
near the water's surface. 

The Naval Undersea Warfare Center (NUWC) Division, Newport, Rhode Island, is 
currently developing an autonomous, 21 inch diameter "torpedo shaped" AUV, known as the 
21UUV, for which near surface operations is envisioned in the future. The 21 UUV shape, in the 
expected 301 inch long version, is shown in figure 1.1. 

The Deep Submergence Laboratory at the Woods Hole Oceanographic Institution. is 
involved in the research and development of the control algorithms for the vehicle, and an 
understanding of the expected environmental disturbances to the AUV will allow a more thorough 
evaluation of the effectiveness of the developed controllers. Also, future decisions concerning 
possible operating regions of the 21] UUV must account for water surface conditions and the effect 


of waves when near surface missions are considered. 
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Figure 1.1 21UUV Profile 


1.2 Research Objectives 


In severe sea conditions, the destabilizing effect of surface waves on an AUV is expected 
to be the limiting factor when considering the upper boundary of useful operating depths 
available to the underwater vehicle. For example, with surface waves in deep water, one would 
expect the amplitude of water motion, and hence the effect of wave forces, to decay with 
increasing depth. Because there is a limit in its ability to stabilize itself, an AUV would have a 
ceiling to its effective operating regime. Therefore, one reasonable measure of an AUV controller 
15 15 performance in the presence of wave disturbances. 

Before the disturbance rejection properties of any controller can be evaluated, the 
properties of the disturbance must be determined. Because the 21 UUV shape is relatively simple, 
existing literature concerning the hydrodynamic forces on similar shaped bodies, i.e. cylinders, is 
abundant. Therefore the first objective of this thesis 1s to apply existing theory to develop a 
model for predicting the forces and moments caused by sea waves on a stationary, slender body 
AUV. Linear wave theory, hydrodynamic strip theory accounting for the precise contour of the 


2] UUV, and the stochastic description of the sea surface are to be used. 





With any model, simplifications of the true physical processes result in model 
inaccuracies. While full scale testing of the yet to be built 21UUV is beyond the scope of this 
thesis, scale model testing in a wave tank is possible. Therefore, the second research objective is 
to conduct tests which either confirm the validity of the wave force and moment model. or 
provide empirical data which allows for the estimation of the hydrodynamic forces and moments 
on the AUV. 

While a precursor of the 21UUV is currently undergoing sea trials which, in part, are 
being used to evaluate controller performance in still water and in steady currents, the testing of 
this vehicle in other than calm sea conditions is a future prospect. Hence, the third research 
objective is to develop a controller for the 21UUV using the same methodology as is expected to 
be used for the actual 21UUV controller, and by simulation, to evaluate the controller's 
performance in the presence of wave disturbances similar to that which might be encountered in 
practice. For simplification purposes, motion in the AUV's longitudinal plane alone is 


considered. 


1.3 Outline of Thesis 


Chapter 2 develops theory allowing for the estimation of wave disturbances on a 
stationary slender body AUV beneath the water's surface. Slender body strip theory and linear 
wave theory are used to develop a method for calculating the transfer function phase and 
magnitude between surface water waves and five of the six forces and moments expected for the 
submerged AUV, both in head and in beam seas. A spectral description of random water waves 
is presented, and using linear time invariant systems theory, a method for calculating the spectra 
of the wave disturbance is shown. Generating a time simulation of waves from their spectral 


representation 1s addressed. 
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Chapter 3 contains a description of the experimental testing performed on a scale model 
of the 21UUV to evaluate the transfer function representation of wave forces and moments 
presented in chapter 2. The experimental apparatus is detailed, as are the series of tests 
performed. Experimental data is compared to theoretical values and largely verifies the earlier 
developed theory. The effect on transfer function magnitude of vehicle forward motion in head 
seas is also investigated. Differences and similarities between the static and dynamic model cases 
are noted. 

Chapter 4 details the general six degree of freedom equations of motion for an 
underwater vehicle. Model simplifications are made accounting for 21UUV body symmetry and 
assumptions concerning maintenance of the vehicle roll angle at 0 degrees. The resulting 
longitudinal plane equations used for subsequent discussion are presented. 

Chapter 5 provides a method of applying sliding control techniques to the 21 UUV in the 
longitudinal plane. Variations of the sliding controller are applied to the vehicle pitch axis, and 
are demonstrated in simulation as an integrated part of the pitch-depth-speed controller. Pitch 
disturbance rejection properties of the controllers is investigated through time simulations, and 
extensions to an adaptive sliding controller are made in an attempt to improve disturbance 
rejection properties. 

Chapter 6 summarizes the results of the thesis and describes the direction of future 


research. 
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Chapter 2 WAVE DISTURBANCE 


The results of linear wave theory provide a first order approximation to the motion of a 
body of water due to surface gravity waves. Some results from the theory are presented here. and 
then are used to develop a model of the stochastic disturbance which could be expected to affect 
an AUV operating near the surface of the ocean, where the wave effects are most prominent. 
This wave disturbance model will be compared to experimental results in a later chapter, and then 
used in AUV dynamic simulations, where the goal will be to reduce the effect of wave 


disturbances through the use of different control schemes. 


2.1 Linear Wave Results 


A more thorough discussion of these results can be found in (Newman 1977), or 


(Faltinsen 1990). 


2.1.1 Regular Waves 


For a single frequency water wave traveling in the direction measured by the angle œQ 
with respect to the Cartesian coordinate frame positive x direction, the free surface elevation 


above the mean free surface can be described by 


C=C, sin(wr — kx cos 0 — kysin Q) (2.1) 
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where G, represents the surface wave amplitude (half the wave height), @ ==* is the circular 
SONA 


frequency, T is the wave period, t is ume, k= is the wave number, and A ts wavelength. 


Wave number, A. is related to circular frequency, €, through the dispersion relation 
q y E p 


= ktanh kh (2.2) 


where g is the acceleration of gravity, and h is the water depth (from mean free surface to ocean 


floor). 


For a wave traveling from one water depth to another, œ remains constant, and the wave 


number, and therefore the wavelength are affected by the change in A. 


Assuming from here on that @ = 0, a water particle's vertical motion, of amplitude €, on 
the surface, decays with depth, z, where z 1s taken positive down, and 1s described by 
ea sinhk(-z+h) sin(Qt — kx) (2.3) 


sinh kA 


while the water particle's horizontal motion is 


= cosh k{-z+h) Des / 
62-6, ——u-cos(mot— kx) (2.4) 


The vertical velocity and acceleraton fields are 


sinhk(—24/ 
w=0€, PEL cos(w1 — kx) er) 
2 sinhk(-z4A) ... n 
a, = -w’ ER sin(wr-kx) (2.6) 


and the horizontal velocity and acceleration fields are 
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shki=-z+h) > 
ዘ= 035, ጸን (637 -- ዴነ) (2:05) 


sinh kA 


coshk(-z+h) a SER = 
= 0 E, HE cos(wr— kx) (2.8) 


The dynamic pressure field in the water column 1s 


cosh k(- coshk(-2+h) Te 
— pg GR an cosh kh sin(w/— kx) (2.9) 


where p ıs the water density and g is the acceleration of gravity. 


2.1.2 Statistical Description of Waves 
In practice, linear wave theory 1s used to simulate irregular seas by the superposition of a 
large number of regular waves (Faltinsen 1990). For a long crested, irregular sea with waves 


traveling in the positive x-direction, the sea surface elevation can be described 


N 
ከ=ኔ 4 sin(Q 1 — k,x € £;) (2210) 


J=1 


where A, œ, k, and e; are the wave amplitude, circular frequency, wave number, and random 


phase of the j-th wave component respectively, and N is the number of wave components used in 
the simulation. The random phase angles are uniformly distributed between Q and 27 radians, 
and the wave number and circular frequency are related through the dispersion relation. Wave 
amplitude is related to the circular frequency through a single-sided wave amplitude spectrum, 


S*(w), and can be calculated from 


A; =25S*(w,)Aw (2.11) 
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Here, Aw is the increment in @ used in the discrete approximation of the spectrum S*(0). In 
implementing equation (2.10), @; 1s chosen randomly and uniformly in the interval @, to @, + ር) 
= @,,, to avoid the repetition of the expression after 277/A@ seconds. It follows that the horizontal 
velocity and acceleration fields, and the vertical velocity and acceleration fields can be simulated 


in the Same manner, as 


= 


coshk (-z+h) . 
= Å . ——— g I-—kk.x 3 17 
ዘ yo A, E sin(Q jt k,x *£,) (22) 


z 


^ coshk,(—z+h) 


= pa ይመ =k 2 

a, (A, nee cos(Q I kx * £,) (2.13) 
yal 
2 inhk,(-2+h) 

Zu sinhk ,(-2+h ው ^ 

wa A en ር05((9 / k;x * £,) (2.14) 
j=l 
N 


= 24 sinhk,(-z+h) . LM ደ 
= 27094 sinhk h sin(@ ፤ ርክ Kx ag E€) (2.15) 


respectively, and that the dynamic pressure can be sımulated by 


Lx TIR h sin(Q t — kx - £;) (2.16) 


The single-sided spectrum S*(@) 1s commonly used in Ocean Engineering applications, 


and 1s defined 


S*(@) = 


che — JWT > 
፡ 15%(፡ dt (020 (2.17) 


0 ዐ<0(0 
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where 
E 
u il 
Bem Jun SED (2.18) 


= ECN EE +) 


is the sea surface elevation autocorrelation function. It can be recognized that S*(0) is related to 


the familiar definition of the power spectrum, D.. (0) (the Fourier transform of the 


autocorrelation function) as 


=D. (0) 020 


S (ys (2.19) 
0 


0)<0 


S*(0) can be calculated in the method described above, that 1s, by first calculating the 
autocorrelation function of a set of wave data and then computing its Fourier transform. The 
assumption made is that sea waves can be described as a stationary random process over some 
short period of time on the order of a few hours. By curve fitting some function of frequency to 
the resulting empirical data, many oceanographers have compactly described the frequency 
content of their data by an empirical formula representing a continuous wave spectrum (St. Denis 
1969). 

The forms of the function used to curve fit wave record data to describe a spectrum are 
various. Bretschneider is credited with proposing the first easily usable two parameter spectrum 
representing seaways in all states of development (Chryssostomidis 1974). The 15th 
International Towing Tank Conference (ITTC) recommended a spectrum of the Bretschneider 


form as the standard international spectrum when information concerning typical sea spectra for 


2 





a specific region of the seas is not available. For seas not limited by fetch, the ITTC 


recommended Bretschneider spectrum has the form 


173 H; 269] 
soe] (2.20) 
o ሀ) 


Here, T, 1s the average wave period, and H, 1s the significant wave height, defined as the average 
n 


of the highest one third of all the waves (15th ITTC 1978). 

The term "sea state” ıs commonly used to describe sea surface conditions ranging from 
glassy seas (sea state 0) to those encountered during hurricane conditions (sea state 9). Using 
data published in (Berteaux 1991) relating sea states to the two parameters above, the ITTC 
recommended spectrum for conditions spanning sea states | through 3 is depicted in figure 2.1. 
It can be seen that as the sea state becomes rougher, the spectrum becomes more peaked, and the 
modal frequency decreases. Also, the majority of the spectrum power is seen to be in frequencies 
below 3 rad/sec, even for the calmest of seas. 

The ITTC recommended spectrum will be selected as the sample wave spectrum in all 
following discussion and simulations. While this spectrum may not be the best available model 
of the actual wave spectrum for a specific application, it ıs assumed to be sufficiently 


representative of the developed model wave spectra for the purpose of this discussion. 
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Figure 2.1 ITTC Spectrum for Seas not Limited by Fetch 
and Conditions Ranging from Sea States 1 to 3 


2.2 Force Predictions 


2.2.1 Load Regimes 
Morison was the first to propose that the horizontal force per unit length on a stationary 


vertical cylinder in waves can be written as 


dF = (on 22 Cua * p2 Czulul) dl (2.21) 


where p is the water density, D is the cylinder diameter, / is the cylinder length, a, and u are the 


horizontal acceleration and velocity of the water at the depth of the cylinder section, and C,, and 
Cp are coefficients which can be determined experimentally (Morison, et al 1950). It is seen that 


this formulation represents two types of forces on a submerged cylinder, the first term 


25 





representing an inertial force proportional to the acceleration of the water at the depth of interest, 
and a second, nonlinear drag term proportional to sign velocity times square velocity of the water 
at the depth of interest. In practice, C, and Cp are dependent on several parameters such as the 
Reynolds and Keulegan Carpenter numbers of the flow, and the surface roughness of the 
cylinder. 

It is therefore possible that in a particular type of flow that either the inertia or drag force 
is predominant. Such is the case for vertical pilings penetrating the water's surface, and it is 
known that the ratios of wavelength and waveheight to cylinder diameter are key parameters in 
predicting the load regime of the waves on the cylinder (Faltinsen 1990). Figure 2.2 depicts these 
load regimes. 

For a Stationary object in a simple harmonic oscillating flow, the time varying total force 


can then be expressed as 


F, = F,sinor|sin@r|+ F, cosor 


ምጫ 
bh 
D 
D 
Nur 


where Fp and F, represent the maxima of the drag and inertia force components, respectively. [t 


can be shown that (Dean and Dalrymple 1984) 


F, 2F, & F, 


F- F (2, 
ü fn pnm 


[-ጋ 
tN 
Uy 


The significance of equation (2.23) 1s that the maximum force on the body 1s not affected 
by additional drag force until the amplitude of the drag is at least one half that of the inertia 
force. For harmonic oscillating flows, such as that caused by regular waves, while even small 


amounts of drag may be important when considering the shape of the load function on a 
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Figure 2.2 Load Regimes ona Vertical Cylinder 
(Adapted from (Faltinsen 1990)) 


Stationary body, the peak amplitude of the force is only affected when the drag component 1s 
greater than one half the inertia force. This implies that if the peak of the regular wave force 1s 
the main concern for a particular submerged body, considering figure 2.2, water particle motion 
with amplitude greater than 0.5 diameters, and perhaps up to 2.5 diameters would produce a peak 
force only as high as the peak force due to the inertia term from equation (2.22). 

The same concepts discussed above will be used to predict the predominant forces on a 
stationary horizontal cylindrical body (the 21 UUV) under waves. When waves cause the motion 
of a water particle at an AUV's depth to be of the order of one UUV diameter or less, it 1s 
expected that the predominant hydrodynamic force on the UUV due to the wave disturbance 
would be inertial in nature. Because an AUV may be deeply submerged, it is not the surface 
wave height to AUV diameter ratio which is of concern, but more appropriately twice the 


amplitude of water particle horizontal or vertical motion at the vehicle depth (which is analogous 





to wave height) compared to the cylinder diameter. For instance, if yaw moment on the body ıs 
of concern, horizontal water motion tangent to the longitudinal (x) axis of the AUV should be 
considered as this 1s the flow which causes the yaw moment. 

While the above analogy is approximate in nature, it provides a means to predict which 
hydrodynamic forces may be of concern when predicting the total load on a cylindrical AUV 
caused by waves. Experimental data will be presented in a later chapter which tests the validity 


of these arguments. 


2.2.2 Inertia Dominated Flow 

Figure 2.3 depicts the axes, force and moment conventions for an AUV used in this and 
subsequent discussions. The body-fixed axes are labeled x, y, and z, with forces X, Y, and Z 
positive in the corresponding positive axis direction, and moments K, M, and N are positive using 
the right hand rule. 

For the long, streamlined body of the 21 UUV, strip theory can be used to calculate the 
hydrodynamic forces and moments imposed on the stationary body by water flow perpendicular 
to the longitudinal axis of the body. Considered below are two body-to-wave orientations, both 
for the horizontal vehicle in an inertia dominated wave force regime. Strip theory is first applied 
to estimate the heave force, Z, the pitch moment, M, and the surge force, X, on the vehicle in 
direct head seas, where wave propagation is perpendicular to the AUV longitudinal axis. Then Y 
(sway force), Z, M, and N (yaw moment) are estimated for the vehicle in direct beam seas, i.e., 
when wave propagation is parallel to the vehicle longitudinal axis. 
2.2.2.1 Head Seas 

To better understand what to expect for Z and M on the AUV body under head seas, first 


considered is a right cylinder of constant diameter equal to the maximum diameter of the 21 UUV, 
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Figure 2.3 Body-Fixed Axis, Force and Moment Conventions for a UUV 


and of the same length as the vehicle. Considering the nearly cylindrical shape of the 21UUV, 
this cylinder model allows for a closed-form solution which roughly approximates the more 
refined solution developed later using numerical methods and taking into account the precise body 
contour of the AUV. 

The vertical force (positive downward) on a stationary horizontal cylinder of length L 
under waves traveling in the negative x direction in an inertia dominated force regime 1s 


calculated here using strip theory as 


Z,(0) » - | Kuss Gc zr) di (2.24) 
L 


where Ka -lnpD'C u. lo determine C,,, the Keulegan-Carpenter parameter is considered and 


15 found for vertical water motion as 


KC, =w, T/D ( 


ho 
n2 
ፍዳ 
— 


2፡7 





where 


EE sınhk(-7+h) ነ 
Wan = 06, sınhkh (2.26) 
is the amplitude of the vertical velocity from equation (2.5), and T is the period of the harmonic 


wave. It is seen that 


KC -2mnC-/D (2.27) 


where €, 15 the maximum vertical displacement of the water particle from its neutral position. It 
is the assumption here that ¢,/D, the "displacement" parameter, is 1 or less and that the resultant 
hydrodynamic force is inertia dominated with C,, = 2 (Dean and Dalrymple 1984). 

Returning to equation (2.24), vertical water particle acceleration is taken at the centerline 


depth of the cylinder. Using linear wave theory and recalling that the wave is now traveling in the 


negative x direction, a (x, z, 1) taken from equation (2.6) can be expressed 


a,(x, 2,1) = A,(z)sin(kx + wr) (2.28) 


Then, recalling that the force Z is taken positive down while the wave elevation 1s taken positive 


up, 
L2 
Zell) A | sin(ke + Wt)dx 
Br (2.29) 
= —K,,,A,2sin &sin(wr) 
and 
Z| =|KysAs2 sine (2.30) 





28 





Similarly, M about the mid-length position on the same stationary horizontal cylinder 


under the same waves 1s calculated using strip theory as 


L2 


M, (t) 2 Kai 4A, J xsinQoc e on)dx a 
-LI2 (2.31) 
Kid uM. | 
= i (2sin + —kLcos*+)cosa@t 
and 
ees eal l 
IM,| esse T COSE. (232) 





The surge force, X, can be estimated by calculatıng the difference ın force between the 


back and front ends of the cylinder due to the difference in the undisturbed dynamic pressure: 


ዊዊ ] 2 coshk(-7+h) = kp 
11011. በርን መን 1 COS) 


oe 2 cosh k(-2+h) aze kL 
Xi ዘ d npgD a cosh kh sin 2 


For a given AUV depth in the water column, the above formulation relates wave number, 
k, to the magnitude of Z, M, and X for unit amplitude surface waves. Since wave frequency is 
directly related to wave number by the dispersion relation (equation (2.2)), equations (2.29), 
(2.31), and (2.33) can be used to solve for the magnitude and phase of the transfer function from 
C, to Z, M, and X. 

For a more refined estimate of Z, M and X, the cylinder model of the 21UUV is 
abandoned and the precise body contour of the AUV is accounted for. Then, from equation 


(2.24), 


L/2 
2 (0=-A, | K y Cosina + 00d (2.34) 


-L/2 
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where እ, (x) =2pxrC Dx). C,, =2 and constant is still assumed, and using numerical 
methods with a look-up table for D(x), the amplitude of Z, can be found for all œ and arbitrary 
phase. Similarly, 

L/2 


My(t)= Ay [x Ky(x)sin (kx + or )dx (2.35) 


=Ur2 
where the same method can be used to find the amplitude for M, for all œ and arbitrary phase. 


Calculating X requires the integration of the dynamic pressure over the vehicle contour at both 


ends, namely 


R R 
X,G) 22n]| pp Gs C) rdr - 2n | Pye) (2.36) 
0 0 


As examples of the calculated transfer function magnitudes, figures 2.4, 2.5, and 2.6 


9፪0፡ን]፤ M, o)|/ E, (0)| and [X (o)|/ 











compare predicted IZ, (ዐ))/ 





።,.(03) for the 21UUV in head 


seas, 30 meter deep water and at various depths using the two methods described above. 

Figures 2.4, 2.5, and 2.6 show decreased transfer function magnitudes with increased 
depth of the cylinder, as would be expected due to the decay of water motion with depth. Also 
observed in these figures is a shift of the peak of the magnitude of the transfer functions to lower 
frequencies with increased depth of the vehicle. This can be explained by realizing that higher 
frequency waves decay more rapidly with increasing depth in the water column than do lower 


frequency waves. 


30 





IZi/Zetal for 21 UUV in 30 meter Water, Head Seas 









2500 
z-5m 
solid - cylinder 
2000 dash - body contour 1 
፳ 1500) 
= 
"e 
D 
N 
N 1000 | 
500 
0 re d 
0 2 35 


wave frequency (rad/s) 


Figure 2.4 Heave Force Transfer Function Magnitude for 21UUV in Head Seas 
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Figure 2.5 Pitch Moment Transfer Function Magnitude for 21 UUV in Head Seas 
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Figure 2.6 Surge Force Transfer Function Magnitude for 21 UUV in Head Seas 


2.2.2.2 Beam Seas 

In left beam seas, the AUV longitudinal axis is considered to be rotated 90% from the 
incoming wave direction, and referring to figure 2.3, the regular wave propagation direction is 
taken in the positive y direction in the body-fixed coordinate system. Considering the range of 
wavelengths over the range of wave frequencies which are of interest, it is noted that A/D 2 13 
for all wave frequencies below 3 rad/sec, and the approximation of uniform water acceleration 
across the diameter of the AUV is made. Strip theory then allows for the calculation of Y, Z, M, 
and N for the AUV in beam seas, while the roll moment, K, though expected to be of 
significance, cannot be reasonably calculated in this manner. Experimental methods best allow 
for determination of K in beam seas, and these will be explored in a later chapter. 

Because Y and N are caused by horizontal water motion, the Keulegan-Carpenter 


number considered is that in the horizontal plane, namely 


97 








KC, > ር) (2.37) 


where 


En cosh k(-2+h) ን 
x no - sinh kh (2.38) 


Here the horizontal displacement parameter &,/D<1 resulting in inertia dominated 
hydrodynamic forces and C,, =2 is assumed. 

Where the previously used cylinder model of the 21UUV can be used to calculate Y and 
Z, using this approach to calculate M and N would predict zero moment about the mid length 
position of the AUV, and therefore only the body contour method is used to calculate M and N in 
beam seas. 


Using the two methods previously described, 


Y (t)= A, K p Lcosœt (2.39) 
using the cylinder model of the body, or 


LI 
Y,(r) » A, coser. | Ky Gd (2.40) 


-L/2 


using the body contour of the vehicle to calculate K,,,(x). Here A, — A, and is taken from 
equation (2.8), Ky. = Ky: and K,,,(x) 2 K,4(x) due to the symmetry of the vehicle. Similarly, 


Z is calculated 


Zg(1) » - A,K y 4Lsin Ot (2.41) 
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or 
Er? 


#,(/)= -4, ጓክ ወያ | ጸ,.()ዐ) (2.42) 


-L/2 


for the cylinder model and body contour model, respectively. 


The moments M and N are found from 


ET 
M(t) = Asinar | x Ky, GOdx (2.43) 


-L/2 
and 


ያጋ 
Na(e)=A, cost |!። ዳ/,/;(፡)ዐ፡ (2.44) 


E 
respectively. 
Figures 2.7 through 2.10 show examples of the calculated transfer function magnitudes 


for |Y,(0)//%, (0), |Z,¿(0)/, (O), [M,(0)/5, (0), and |V,(@)|/|C,(@)|, respectively for the 

















21UUV in beam seas, 30 meter deep water and at various depths using the methods described 
above. Comparing figures 2.7 and 2.8, the transfer function magnitudes are identical for the Y 
and Z forces except at low frequencies where, because of the larger wavelength to water depth 
ratio, the water particle motion decays more rapidly with depth for the vertical motion than for 
horizontal water particle motion. The same can be said when comparing the M and N moments 
in figures 2.9 and 2.10. Because there is little fore-aft asymmetry in the 21UUV, the predicted 
pitch and yaw moments in beam seas are seen to be relatively small compared to the predicted 
pitch moment in head seas (figure 2.5). Finally, comparing Z,, in head seas versus Zp in beam 


seas (figures 2.4 and 2.8), it is observed that at low frequencies, the magnitudes of the transfer 
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i 
l 
| 


mum 


functions are identical. At higher frequencies where the wavelength 1s shorter and of the order of 
the vehicle length, [Zala is predictably smaller than 1.1.1 

The application of strip theory in calculating forces and moments on an AUV in head and 
beam seas in an inertia dominated hydrodynamic force regime has allowed for the prediction of 
the transfer function from surface wave amplitude to forces and moments on the AUV. The 
standing assumption has been that water particle motion at the depth of the AUV is small enough 
so that nonlinear hydrodynamic form drag 1s insignificant when compared to the linear 


hydrodynamic inertia forces. 
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Figure 2.7 Sway Force Transfer Function Magnitude for 21 UUV in Beam Seas 
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Figure 2.8 Heave Force Transfer Function Magnitude for 21UUV in Beam Seas 
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Figure 2.9 Pitch Moment Transfer Function Magnitude for 21 UUV in Beam Seas 
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Figure 2.10 Yaw Moment Transfer Function Magnitude for 21 UUV in Beam Seas 
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2.2.3 LTI Systems with Stochastic Inputs 


Figure 2.11 depicts a linear time invariant (LTI) system with stable, proper transfer 


function G(s), input u(s), and output y(s), where s = (4) is the Laplace operator. 


u(s) G(s) y(s) 


Figure 2.11 LTI System 


Such systems have long been studied, and presented below is a well-known result which will be 
used in later discussions. For a thorough treatment of the subject of LTI systems with stochastic 
inputs, the reader can consult (Papoulis 1984), wherein the proof of the following result is 
contained. 

For the system in figure 2.11, if u(t) is a known input, and G(t), the impulse response of 


the transfer function G(s) is known, then y(t) is known and can be expressed 


y(t) = y(0)+ [Ge 2) dt (2.45) 


0 


Considering the case when u(t) is a stationary, random process with known power 
spectrum, then y(t) is also a stationary, random process. The power spectrum of y(t) can be 


calculated as 


o. (0) 2» |GCjo)| $,, (0) (2.46) 


uu 
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where IGCjo)| is the magnitude of the transfer function G(s) evaluated at s = j@. It follows then, 


that the single-sided spectrum of y(t), from equation (2.19) is 


" lo. (0) 20 
s (0), 92 (2.47) 
0 O<0 


The significance of the above result is that the previous discussion relating surface wave 
action to forces and moments on an AUV has been cast in such a framework. It can be seen that 
if the magnitude of the transfer functions between sea surface waves and forces and moments on 
an AUV are known, then the statistics of the forces and moments on the AUV body can be 
determined. 

As an example, the pitch disturbance spectrum on the 21UUV in head seas and sea state 
2 conditions are calculated for various depths of the vehicle and depicted in figure 2.12. In 
generating these spectra, the ITTC recommended Bretschneider wave amplitude spectrum for sea 
state 2 was used, as well as the body contour generated transfer function from wave amplitude to 
pitch disturbance depicted in figure 2.5. 

Figure 2.13 depicts two possible time realizations of this pitch disturbance which are 
generated using the technique described in section 2.1.2 for generating time realizations of 
surface waves. The differences between the two realizations are due to the random phase used in 
each simulation. There are, of course, an infinite number of possible realizations of this pitch 


disturbance, each having the spectral representation depicted in figure 2.12. 
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Figure 2.12 Pitch Disturbance Spectrum for 21UUV in Head Seas 
for Sea State 2 and at Various Depths 
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Figure 2.13 Two Possible Time Realizations of Pitch Disturbance 
to 21UUV in Head Seas, Sea State 2 and at 10 Meters Depth 
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Chapter3 EXPERIMENTAL TESTING 


In this chapter, experiments which test the theory of chapter 2 are discussed, and results 
of tests conducted on a 21UUV model are presented and compared with the earlier developed 
theory. While chapter 2 theory deals with forces on a stationary body, the wave forces on a 
forward moving AUV are also of interest as many AUV missions are conducted while the vehicle 
is moving with forward velocity. Also presented here, then, are experimental results of wave 


forces on a forward moving AUV model. 


3.1 Experimental Setup 


The experimental apparatus and AUV model are shown ın figures 3.1 and 3.2. 


3.1.1 Experimental Apparatus 

The Massachusetts Institute of Technology's Ocean Engineering Testing Tank was used 
to conduct model testing. The tank has dimensions 110 feet (length) by $ feet (width) by 4 feet 
(depth), is filled with fresh water, and is equipped with a wave maker and moving carriage. 

The carriage assembly is suspended by rollers from a cylindrical beam fixed to the 
ceiling along the length of the tank. The carriage, on which a mast and AUV model were 
mounted, are capable of sliding the length of the tank, with the AUV model submerged in the 
tank water. Also affixed to the carriage assembly is a belt drive which can propel the carriage at 
speeds up to 2 meters per second along the beam. The speed of the carriage ıs controlled, and 


from a dead stop, the belt drive reaches a desired speed within 2 seconds of activation. 
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Figure 3.1 Experimental Apparatus Setup 
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Figure 3.2 21UUV Model Used in Testing 
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The wave maker, located near one end of the tank, consists of a rigid metal wall spanning 
the width and depth of the tank. The metal wall is allowed to pivot about its attachment to the 
bottom of the tank, and it 1s driven by a hydraulic actuator mounted at its top. Waves of 
frequencies between 0.2 Hz and 3.0 Hz can be generated. At the far end of the tank from the 
wave maker is densely packed plastic netting suspended in the water which acts as a wave 
suppresser. The suppresser absorbs much of the wave energy as it reaches the "beach" end of 
the tank, thus largely reducing the amount of reflected wave energy in the tank. 

The wave probe used for measuring wave height uses two parallel copper wires 
separated by approximately one centimeter mounted on a suff frame and positioned vertically in 
the water. A potential 1s applied between the two wires, and the varying resistance, resulting 
from the change in water level due to waves, 1s the means by which water elevation is measured. 


The wave probe was calibrated at the beginning and end of each data collection set. 


3.1.2 AUV Model 

The AUV test model was manufactured as a 1:4.188 scale model of the 301 inch long 
version of the 2] UUV being developed at NUWC. The model body contour is precisely that of 
the 21UUV, including the contour of the tail section and fins. The model was constructed in 5 
parts, and then assembled. The nose and tail sections were manufactured from PVC, while the 
two inner cylindrical sections were made from hollow cast acrylic tubing. The sensor section of 
the model, manufactured from 6061-t6 aluminum, housed a 6-axis strain gauge sensor which was 
mounted to the model at the sensor's bottom and to the rigid support mast at the sensor's top. As 
a result, the resultant hydrodynamic forces and moments on the AUV model were transmitted 
through the sensor to the rigid support mast, allowing for their measurement. The five sections 
assembled as depicted in figure 3.2 and resulted in a streamlined model of the full scale 21 UUV. 


The 6 axis strain gauge sensor was calibrated the first and last days of model testing. 
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During data collection, the sensor's 6 channels and the wave probe's 1 channel were 


simultaneously sampled at 30 hertz, with the data being recorded by a 386 personal computer. 


3.2 Testing 


3.2.1 Overview 
Three series of tests were conducted: two series where the AUV model was kept 
stationary, and the third where the model was towed through the water with forward speed. In 
the first group of tests, the stationary model was oriented with its longitudinal axis perpendicular 
to the oncoming wave crests, 1.e., as if in head seas. In the second series of tests, right beam seas 
were investigated and the stationary model was oriented with its longitudinal axis parallel to the 
oncoming wave crests. In the third series of tests, the model was towed with a fixed forward 
velocity counter to the direction of the wave propagation, simulating an AUV underway in head 
seas. For the tests involving a stationary model, the wave gauge was positioned to measure the 
water elevation at the mid-length position of the model, thus allowing for phase comparisons 
between the wave elevation and the forces and moments on the model. 
The parameters varied during the course of the testing were: 
(1) wave amplitude, 
(2) wave frequency, 
(3) AUV speed and orientation, and 


(4) AUV depth 


3.2.2 Scaling Considerations 


3.2.2.1 Wave Frequencies 
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It is shown in chapter 2 that the ratio of wavelength to AUV length is a primary factor 


in the transfer function between the surface wave motion and the forces which affect the AUV. 


In addition, figure 2.1 depicts the range of wave frequencies over which the majority of wave 


energy is expected for a variety of sea conditions. Therefore, the frequencies of waves generated 


during the tests were chosen such that they produced a wavelength similar in scale to the mode! 


21UUV length as full scale waves would produce relative to the full scale 21UUV length in a 


similarly scaled water depth. 


An example clarifies the calculation: 


Example of Wave Frequency Scaling 

Given: 
Wave tank depth 
Scale of model 
Full scale wave frequency (for example) 
Full scale 21UUV length 
Gravity 

Calculation: 
Full scale depth (model depth / scale) 
Full scale wavelength (equation (2.2)) 
Model wavelength (scaled) 


Model frequency (equation (2.2)) 


48 in (1.2192 m) 
1:4.188 

2 rad/sec 

301 ın (7.6454 m) 


9.806 m/s? 


5.106 m 
14.984 m 
3.578 m 


4.093 rad/sec 


Table 3.1 contains the frequencies and wavelengths of waves (full scale and resulting 


model) used during testing. The 20 frequencies of full scale waves indicated in table 3.1 span the 
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range of frequencies expected of ocean waves as described by the ITTC recommended wave 
spectrum. 
3.2.2.2 AUV Speeds 

While conducting the tests during which the model AUV was towed, the Froude number 
of the full scale 21UUV was considered in determining the velocity at which to tow the model. 


Froude number similitude implies 








O, 
በ L = Im = Fr, 
gL, p (3. | ) 
= T 
U pU U, 24, 


where fs and m represent full scale and model, respectively. Data was collected at the two model 
tow speeds shown in table 3.2, and while higher tow speeds were considered, sensor load capacity 


precluded higher speed testing. 


46 



























Table 3.1 Frequencies and Wavelengths Investigated During Testing 
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Table 3.2 Tow Speeds Investigated During Testing 
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3.2.3 Tests Conducted 


Table 3.3 summarizes the 320 trials conducted during the course of testing. 


Model Wave Amp # of 
Centerline per Frequencies 
Depth (m) Frequency 





Table 3.3 Summary of Tests Conducted 


The 20 wave frequencies referred to in table 3.3 are those listed in table 3.1. 


3.3 Test Results 


3.3.1 Raw Data 

The seven channels simultaneously recorded during each of the trials included the six 
axes from the sensor mounted inside the model body plus the wave gauge output. An example of 
the seven channels sampled (with force, torque, and wave amplitude conversions applied) during 
one test run is depicted in figures 3.3a through 3.3d. This particular sample produced three data 


points for the case of beam sea waves of 5.321 rad/s for the 0.379 meter deep model. As ıs 
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shown in figure 3.3, three wave amplitudes were generated during each trial when the model was 
held stationary. During tests in which the model was towed, wave amplitude was held constant 


during the course of each data run. 
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Figure 3.3a Sample of Surge Force and Roll Moment Raw Data. Beam sea effects are 
investigated here, and during this data collection run, wave amplitude was increased in three 
distinct steps as shown in figure 3.3d. While the amplitude of surge force, X, 1s only slightly 
larger than the sensor and A/D converter resolution, roll moment, K, more fully spans the sensor 


and A/D converter full range. 
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Figure 3.3c Sample of Heave Force and Yaw Moment Raw Data 
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Figure 3.3d Sample of Water Elevation Raw Data 


3.3.2 Signal Processing 


The frequency of encounter between the model and waves ıs given in (Newman 1977) as 


W, = (0$ — KU cosO (3.2) 


where (0, 1s the wave frequency, U is the forward speed of the model, and 0 is the angle between 
the model x axis and the direction of wave travel. The highest frequency of encounter between 
the model and waves during testing was evaluated as 11.0 rad/s, or 1.75 Hz. 

Prior to evaluating the amplitude of the signal coming from each of the seven channels, 
data from each channel was digitally filtered using a Chebyshev type II lowpass, stopband ripple 
filter (MATLAB 1992). The 9 pole filter had a cutoff frequency of 3.5 Hz and a stopband of 
negative 60 dB. A Bode plot of the filter frequency response is depicted in figure 3.4. The 
signals were first filtered in the forward, and then reverse directions to yield a zero phase shifted, 


filtered version of the output signals. 
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Figure 3.4 Frequency Response of Filter 


Two examples of sampled and filtered signals (superimposed) are depicted in figure 3.5. 
The first of the two signals shown is the 35 to 45 second window of the Y force depicted in figure 
3.3b. The second of the two signals shown is a 25 second window of the X force recorded while 
the model was being towed at 0.489 m/s under 2.251 rad/s waves. Filtering a low frequency and 
relatively noiseless signal such as Y in figure 3.5 leaves it virtually unchanged. The signal 
representing X in figure 3.5 has a significant level of high frequency carriage rumble noise 
superimposed upon it, and filtering a noisy signal such as this allowed better estimation of the 


amplitude of the wave induced hydrodynamic force. 
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Figure 3.5 Two Samples of Unfiltered and Filtered Signals 


3.3.3 Experimental Data vs. Theory 

After the signals were filtered, amplitudes of the signals were determined and the ratios 
of force and torque to wave amplitude were calculated and plotted versus frequency. 
Additionally, the phase difference between the wave elevation sinusoid and force / torque signals 
were measured and plotted for the cases when the model was stationary. The vertical and 
horizontal displacement parameters were calculated for each test conducted and were found to be 
less than 0.4 in all cases, establishing the tests within the range of displacement parameters 
assumed in chapter 2. 
3.3.3.1 Stationary Model in Head Seas 

Figures 3.6 through 3.8 compare theoretical and experimental transfer function 


magnitude and phase information for X, Z, and M versus wave frequency for the model at a 
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centerline depth of 0.379 m. The theoretical curves were developed using the methods described 
in chapter 2 with the body contour of the model taken into account. Similarly, figures 3.9 
through 3.11 compare the same transfer function magnitude and phase information for the model 
at a centerline depth of 0.787 m. 

Resonance in the wave tank across its width at wave frequencies of 3.3 and 5 rad/s 
appear to cause erratic data near these frequencies during the course of testing, and the result is 
seen in the data presented here. 

While the shallow and deep model data presented for the X and M transfer functions is 
well predicted by the theory both in phase and magnitude, the method used to predict the transfer 
function for Z fails to include a force component which accounts for the resultant vertical force 
when the water wavelength is the length of the model body. The predicted zero in X is observed 
at or near this frequency, as both figures 3.6 and 3.9 show in the magnitude and phase plots The 
phase of this unpredicted Z force is consistent with that which would be expected of vertical drag 
proportional to wave velocity in the aft section of the model body. The inclusion of such a drag 
component into the heave force model was investigated, and produced a far worse low frequency 
fit to the data than that presented in figures 3.7a and 3.10a. Because model accuracy is deemed 
more important at lower frequencies where the majority of wave spectral energy is expected, the 
previously developed model for heave force will be used in subsequent discussion. 

It is seen that X, M and Z are all approximately linearly related to the wave amplitude 


C(@), thus allowing for their transfer function representation. 
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Figure 3.6a Surge Force Transfer Function Magnitude for Shallow Model in Head Seas 


Note: The data points 'x', 'o', and 'x' represent data taken from 
the model under waves of increasingly higher amplitudes, respectively. 
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Figure 3.6b Surge Force Transfer Function Phase for Shallow Model in Head Seas 
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Figure 3.7a Heave Force Transfer Function Magnitude for Shallow Model in Head Seas. The 


erratic effect of wave tank resonance across its width at 5 rad/s 1s clearly seen. 
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Figure 3.7b Heave Force Transfer Function Phase for Shallow Model in Head Seas 


56 





IMI/IZeta! (N-m/m) 


IMl/IZetal for Stauonary Model in Head Seas 


X 


180 Z= 0.379 m = 
solid - theory 
160 points - data 


140 | o 
120 a 5 
x 

100 D 

80 5 

60 iS x 

x 
40 e x 


20 ኔ 


2 3 4 sa MEC 7 E 


wave frequency (rad/s) 


Figure 3.8a Pitch Moment Transfer Function Magnitude for Shallow Model in Head Seas. 


Again, the effect of wave tank resonance at 5 rad/sec is seen. 
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Figure 3.8b Pitch Moment Transfer Function Phase for Shallow Model in Head Seas 
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Figure 3.9a Surge Force Transfer Function Magnitude for Deep Model in Head Seas 
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Figure 3.9b Surge Force Transfer Function Phase for Deep Model in Head Seas 
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Figure 3.10a Heave Force Transfer Function Magnitude for Deep Model in Head Seas 
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Figure 3.10b Heave Force Transfer Function Phase for Deep Model in Head Seas 
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Figure 3.11a Pitch Moment Transfer Function Magnitude for Deep Model in Head Seas 


Phase of M Relative to Surface Wave, Head Seas 





200 
| 2= 0.787 
I50r solid - theory 
points - data 
100 - 
50 
=ዕ 
ረጋ 
= 
6) 0 
d 
A 
መጃ 
-50 , 
0 a ሽ . x a 
-100 9 p A A x e » a e i 
x X * e X 
-150 
-200 
1 2 3 4 5 6 7 8 


frequency (rad/s) 


Figure 3.11b Pitch Moment Transfer Function Phase for Deep Model in Head Seas 
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3.3.3.2 Stationary Model in Beam Seas 

Figures 3.12 through 3.16 compare theoretical and experimental transfer function 
magnitude and phase information for Y, Z, K, M and N versus wave frequency for the model at a 
centerline depth of 0.379 m. The theoretical curves were developed using the methods described 
in chapter 2 with the body contour of the model taken into account. Similarly, figures 3.17 
through 3.21 compare the same transfer function magnitude and phase information for the model 
at a centerline depth of 0.787 m. 

Again the data 1s well predicted by theory in phase and magnitude, with the exception of 
the prediction of K, where no theory is presented. The scatter in data for M and N is noted, as 
are the relatively small values of M and N when compared to M in head seas. Because the 
magnitudes of the M and N moments in beam seas are relatively small, sensor axis crosstalk and 
the effect of imperfect waves are possible causes of the scatter in data depicted in figures 3.15, 


ESO 5:20: and 3.21. 
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Figure 3.12a Sway Force Transfer Function Magnitude for Shallow Model in Beam Seas 
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Figure 3.12b Sway Force Transfer Function Phase for Shallow Model in Beam Seas 
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Figure 3.13a Heave Force Transfer Function Magnitude for Shallow Model in Beam Seas 
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Figure 3.13b Heave Force Transfer Function Phase for Shallow Model in Beam Seas 
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Figure 3.14a Roll Moment Transfer Function Magnitude for Shallow Model in Beam Seas 
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Figure 3.14b Roll Moment Transfer Function Phase for Shallow Model in Beam Seas 
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Figure 3.15a Pitch Moment Transfer Function Magnitude for Shallow Model in Beam Seas 
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Figure 3.15b Pitch Moment Transfer Function Phase for Shallow Model in Beam Seas 
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Figure 3.16a Yaw Moment Transfer Function Magnitude for Shallow Model in Beam Seas 
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Figure 3.16b Yaw Moment Transfer Function Phase for Shallow Model in Beam Seas 
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Figure 3.17a Sway Force Transfer Function Magnitude for Deep Model in Beam Seas 
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Figure 3.17b Sway Force Transfer Function Phase for Deep Model in Beam Seas 
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Figure 3.18a Heave Force Transfer Function Magnitude for Deep Model in Beam Seas 
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Figure 3.18b Heave Force Transfer Function Phase for Deep Model in Beam Seas 
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Figure 3.19a Roll Moment Transfer Function Magnitude for Deep Model in Beam Seas 
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Figure 3.19b Roll Moment Transfer Function Phase for Deep Model in Beam Seas 
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Figure 3.20a Pitch Moment Transfer Function Magnitude for Deep Model in Beam Seas 
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Figure 3.20b Pitch Moment Transfer Function Phase for Deep Model ın Beam Seas 
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Figure 3.2la Yaw Moment Transfer Function Magnitude for Deep Model ın Beam Seas 
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Figure 3.21b Yaw Moment Transfer Function Phase for Deep Model in Beam Seas 
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3.3.3.3 Forward Moving Model in Head Seas 

Figures 3.22 through 3.24 show experimental data depicting transfer function magnitude 
information for the dynamic components of X, Z and M versus wave frequency for the model at a 
centerline depth of 0.379 m and 0.489 m/s forward speed. Theoretical curves for the stationary 
model are plotted with the data for reference. Similarly, figures 3.25 through 3.27 show 
experimental data for the model at the same depth with a forward speed of 0.733 m/s, with the 
same theoretical stationary model curves plotted for reference. 

Continuous shedding of vortices behind the model 1s now expected due to the model 
forward speed. The resulting form drag due to this vortex shedding makes the stationary model 
theory less applicable. Deviation from the results of the static model testing is clearly seen for X 
and Z at higher frequencies, and particularly where the zero of the X transfer function was seen 
at approximately 6 rad/s for the static model AUV. The similarity of the M transfer function 


magnitude for the static and dynamic AUV cases is noted. 
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Figure 3.22 Surge Force Transfer Function Magnitude for Model at 0.489 m/s in Head Seas 
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Figure 3.23 Heave Force Transfer Function Magnitude for Model at 0.489 m/s in Head Seas 
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Figure 3.24 Pitch Moment Transfer Function Magnitude for Model at 0.489 m/s in Head Seas 
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Figure 3.25 Surge Force Transfer Function Magnitude for Model at 0.733 m/s ın Head Seas 
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Figure 3.26 Heave Force Transfer Function Magnitude for Model at 0.733 m/s in Head Seas 
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Chapter 4 AUV DYNAMICS 


4.1 Equations of Motion 


4.1.1 Coordinate Systems 

The motion of an AUV in its body reference frame is typically related to its motion in an 
earth-fixed reference frame versus that of a true inertial frame. For the relatively slow moving 
ocean vehicle, the neglected motion of the earth's surface is of small consequence. 

Figure 4.] depicts the three ordered rotations, known as Euler angles, which describe the 
coordinate transformation from the earth-fixed to body-fixed reference frames. To understand 
the rotations which describe this change in frames, it can be imagined that a body first begins in a 
neutral earth fixed attitude: facing true north, level with the horizon both to the north and east. 
To reach the final vehicle attitude, the first rotation is taken about the normal (z) axis, is called 
yaw and denoted y, and can be thought of as a heading change (positive right). The next rotation 
15 pitch, taken about the transverse (y) axis, and is denoted 0 (positive nose up). The third and 
final rotation is roll, taken about the longitudinal (x) axis, and is denoted © (positive right). The 
order of these rotations 1s not arbitrary, and taken in a different sequence, could result in a 
different set of angles for the same final body attitude. All body attitudes, except those including 
0 2 x90*, can be uniquely described by the rotations of the three Euler angles. The singularity 


associated with 9 = 290" can be avoided by using a four parameter method to 
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Top View 
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Figure 4.1 The Euler Angles 


7] 





describe the three basic rotations, but the natural insight into the body attitude is lost in this 
approach. For this discussion, 16] < 90 is a standing assumption. 
Based on the SNAME (1950) notation, the motion of an AUV can then be described with 


the following vectors: 


n=|n ml O espe 
v= [vi v | vi=[u v w] SE vj (4.1) 
E | Doe Rem E N) 


where T] represents position and attitude information which fixes the vehicle in the earth reference 
frame, v represents translation and rotation information relative to the body fixed reference 
frame, and 1 represents external forces and moments acting on the AUV in the body reference 
frame. 

The scalars x, y, and z from equations (4.1) represent the vehicle position in an earth 
fixed, right hand, three dimensional reference frame with the x-axis and y-axis pointing to the 
horizon, and the z-axis positive downward. The scalars 9, 9, and y are the previously mentioned 
Euler angles of roll, pitch, and yaw, respectively. The scalars u, v, and w represent body-fixed 
reference velocity along the longitudinal axis (surge), the transverse axis (sway), and the normal 
axis (heave). The scalars p, q, and r refer to right hand rule angular rotation rate about the body 
longitudinal (x) axis, transverse (y) axis, and normal (z) axis, respectively. The values X, Y, and 
Z represent forces along the body x, y and z axes respectively. K, M, and N are moments along 
the axes which would cause positive roll, pitch and yaw, respectively. Table 4.1 summarizes the 
SNAME notation. 


The description of the vehicle's flight path in the earth fixed frame given the body fixed 


motion is governed by the transformation matrix 1) (ኻ. ) , 1.6. 
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MOTION FORCES LINEAR POSITIONS 
and and and 
MOMENTS | ANGULAR EULER 
VELOCITY ANGLES 
(n) 


translation in x (surge) 





Table 4.1. SNAME Notation Used for Ocean Vehicles 


Tl - Ji(mv, 
cycO | —sycoó-ccysOso. — syso-c cy cosO 
(4.2) 
Ji(m) 2| sy c8 CV có t+ sosOsy =c sO + sOsy có 


-s0 ርፀናዕ ርፀርዕ 


where s(e), c(e), and 1(e) represent sin(e), cos(e), and tan(e), respectively. 


Similarly, the description of the vehicle's earth fixed rotation given the body fixed 


rotation vector is governed by the transformation matrix J,(n,), 1.e. 


mn. =J,(n,)v, 
| $010 core 
J,(n,)=/0 ርዕ “#0 
0. 40/ ር90- ርዕ / ርፀ 


(4.3) 


uo 





Since each of the two reference frames uniquely maps to the other, it follows that the two 


transformation matrices are invertible, and the reverse transformatıons also hold, ı.e. 


v, 7 J; (ni (4.4) 
v,2J4 (ኻቃ)ኻ; (4.5) 

are also true. 

More generally, it can be written 
ኻ >= 7(ኻ:)ነ (4.6) 
v=J "(nn (4.7) 
where 
7 ¡(n2) 0 | 
J(N,) = (4.8) 
0 3-(ኻ.) 


4.1.2 Rigid Body Dynamics 

The rigid body equations of motion for an underwater vehicle can be written about an 
arbitrary origin. When considering an AUV with body symmetry, it is convenient to take the 
origin as the intersection of the longitudinal, lateral, and normal axes of the body. If the cross 
moments of inertia about the center of gravity are negligible, then using the parallel axis theorem, 


the rigid body equations of motion can be expressed as: 
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mu — vr * wq- ን ከ). ፡5(79 ያ) + 2ራ(#ዐ" + ያ)]= ፳ 


m|v — wp e + DE za (qr — p) - xg(qp* r))- Y 


m[w — uq + vp — zc (p^ tq )txgs(rp-q) * ys(rq + p)=Z 
Lp+(,—1,qr+ ml yo (w—uq +vp)-2¿(v—-wp +ur)]= K 


(4.9) 


1,9+(U,—1Jrp+mlz(u—=vr+wq)-x¿(w-uq + vp)]= M 


Lr, -H)pq m[xg(v — wpetur) - yg(u vr t wq)] 2 N 


where r, = bee Ye x] represents the displacement from the origin to the center of gravity. 


The above equations can be written in the more compact form (Sagatun and Fossen, 1991): 


m 


eo re 


where 


<> 


Mo 


MeV + Cog VIV = Tag 


0 0 0 MZ 
m 0 —MZG 0 
0 m myg ME 
—M2G myg l 0 
0 —MXG 0) i 
"nx. 0 0 0 


(4.10) 


—nmyg 


mx 


(4.11) 


and the Coriolis matrix, though not unique, can be written in a skew symmetric form 
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0 0 0 m ycd t zF) AMOS) 
0 0 0 —M(YGP EW) M(Z¿F+X¿P) 
0 0 0 MEP >V) M3 9 +4) 
Ca(v)= | 
-m(ycQ+Zgr) M(yg p * w) Mtz cp Y) 0 Ir 
m(x.9-w) —IM(2GF + xsp) m(26q +u) EE 0 
#፲[ ጳዕ7 * v) m Yor = u) —M(X¿P + Yc9) lq UD 
(4.12) 


MUA nd | 
>= ያ የነ NU) | 
M(XGP + yaq) 
25g 
[p 
0 


Here it is worthwhile pointing out that the above relationships are Newton's equations for 


six degrees of freedom, and that the vector Tp, =[X Y Z K M NJ represents all 


external forces on the vehicle, such as hydrodynamic added mass and damping terms, restoring 


forces and moments caused by gravity and buoyancy, and any controls and disturbances. 


4.1.3 Hydrodynamic Forces 


While the rigid body equations of motion can be exactly derived for a given body, the 


external force vector, Tpp, is usually obtained by using a combination of theoretical and 


experimental methods, which is to say that an exact solution for Tp, is not available. Fossen 


(1994) uses the notation 
Tor = Tp + Ty + TE + TO 


where 


Tr represents the radiation induced forces, and ıs the vector sum of the 


hydrodynamic added mass, potential damping and restoring forces, 


Ty is the vector of viscous damping forces 
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(4.13) 





ፒራ represents the environmental forces, including those caused by currents and 
waves, and 


tc is the vector of control forces, such as those from thrusters and control 


surfaces. 


Then it can be shown that the first term in equation (4.13) can be represented as 


tg =-M,v-C,(v)v -D,(v)v — -gem (4.14) 


ዥ ሙምኬኬ ሜጥማ ም a ne 
added mass and Conolis potential damping restoring forces 


Here, the matrices which comprise the added mass and Coriolis terms are represented as 


X, X, Ar X; X, X; 
DEL EIS NS: 
Z, Z Z; Z; Z; Z; 
M, =- (4.15) 
K, K; K, K; K, k. 
M, M, M; M, M, M, 
N; N. N. N, N, N; 


and though not unique, C,(v) can be written in the skew symmetric form 
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0 0 0 0 ES 


0 0 0 C 0 MEC 
0 0 0 Bere 0 

C v= 8 F 4 (4.16) 
Í IM O CE 
=C} 0 -C -CY 0. 7.5 
Las IC” 0 zer -cC 0 


where 
CH =-X,u-Y,v-Z,w-Z,p-Z,q-Zr 


CP = Xu + Y v+Y,w+ Y,p + Y,q + Y,r 

CO =X u+ Y, v+Z w+ Z,p+ 2,gt+Z;r 
C$ =-X 4 — X v— X w-— AP Ad A 
C, 2 —-X,u - Y,» -Y,w- Y,;p—-Y;q—- Yr 
CA = Xu + Xy + X,w+ X,p+ X,g + X;r 
Cr =-Xu-Yv-Zw-K.p-M;g-N;r 
C, 2 Xu Y,v t Z,w+ K,p+ M,g+ M;r 


En SAME EV EDER DIR Gr 


A 


The potential damping term from equation (4.14), D,(v)v, results from the generation of waves 
on the free surface due to the motion of the AUV, and can therefore be considered negligible 
when the AUV is sufficiently deep to preclude its generation of waves (Newman 1977). The 
restoring forces are those caused by the forces of gravity and buoyancy upon the vehicle, and 


recalling that the z axis 1s taken positive downward, can be represented 
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(W — B)s50 
—(W — B)c0s6 


-(W — B)c8có 


5 (4.17) 
ፎ(ኻ -(157/ — ypB)cOch +(2¿W - 2,8)ር650 


I ፡፡)፡ ጋ ምህ (1) 
| -(x59W - x4 B)c0s$ — (y W — y¿B)50 


where W and B represent the vehicle weight and buoyancy forces, respectively, and 
r; = be Ya zn is the displacement from the origin to the center of buoyancy. 
The second term of equation (4.13), representing viscous damping, can be written in the 


form 


Ty = -D,(v)v (4.18) 


1. 
viscous damping 


where the matrix D, (v) can be calculated with varying degrees of accuracy by either direct 
experimentation or the use of published experimental data such as that in (Hoerner 1965, 1975) 
for similar shaped bodies. 

In this discussion, the characterization of the environmental force term, T; from equation 
(4.13), will be restricted to the forces caused by surface waves, and is the subject of chapters 2 
and 3. The final term in equation (4.13), the control vector tc, is dependent upon the mix of 


thrusters and / or control surfaces which are available for a specific AUV. 


4.1.4 Body Symmetry Considerations 
In the above discussion, equations (4.10) and (4.13) taken together provide the general 


framework for developing the nonlinear, highly coupled equations of motion for an AUV, with 
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the only assumption thus far being that the cross moments of inertia about the center of gravity 
are considered negligible. Further exploitation of the symmetry of an AUV body can result in 
simplified representations of the rigid body mass and Coriolis matrices from equations (4.11) and 
(4.12), as well as those matrices comprising the hydrodynamic added mass and Coriolis forces, 
restoring forces, and viscous damping terms. 

The top-bottom and left-right symmetry of the 21UUV allows for simplification of the 


hydrodynamic added mass and Coriolis matrices, namely 


A c) 0 0 0- 0 
| 6!) 0 | ጢስን 
COTE MEO Zar 
M == (4.19) 
0-— 0 ሰ 0 
Qu 0 ML. 0-M. 0 
VE LU 0 OV RY. 
and the same representation for C ,(v) as in equation (4.16), where now 
15 16 24 
C, =-Z,w-Z.q C,=Yv+Yr C, —-ZwtZgq 
(EMT mpg TC Exu (4.20) 


Ci =-Y;v-N;r Ci =Zw+Mgq Ci =-K;p 


(Fossen 1994). 


Considering only first order viscous damping terms, the form of D,(v) for this 


discussion is chosen as 
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X le 0 0 0 0 0 
0 Y O0 0 Y 
oe 7, 
D,(v)=- (4.21) 
የ Oe ሰ 
0. 0 M, M, 0 
0 . 0 0O N, 


Also assumed for the remainder of this discussion ıs that for the 21UUV, y¿= 0, which 


further simplifies Cop V), Mp, and g(n). 


4.2 Model Simplifications 


Healy and Marco (1992) suggest that the 6 degree of freedom equations of motion for an 
AUV can be divided into three non-interacting (or lightly-interacting) sets of equations for control 


of speed, steering and diving, each involving the state variables: 


(1) Speed system state: u(t) 
(2) Steering system states: v(t), r(t), y(t), x(t) and y(t) 
(3) Diving system states: w(t), q(t), O(t), and z(t) 


The rolling mode (p(t) and @(t)) is left passive in this approach for their vehicle, the Naval 
Postgraduate School AUVII. Their motivation for this approximation is the limited number of 
actuators on their vehicle, consisting of forward thrusters, stern planes and rudders. 

An approach similar to this is currently being used in the development of a sliding 


controller for a precursor to the 21UUV, the NUWC LDUUV, though the rolling mode is actively 
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controlled for this vehicle. With this controller design, inner loop pitch and yaw controllers are 
used to implement desired attitude control action calculated in the outer loop depth and heading 
controllers, respectively. With this model simplification and accounting for the body symmetry 
of the LDUUV which is similar in shape to the 21UUV, control of the vehicle in the dive and 
steering planes are seen as nearly identical problems, the only difference being the restoring 
forces which act in the dive plane and are absent in the steering plane. The speed system and roll 
system controllers are somewhat more straightforward to design and implement because of the 
actuators available on the vehicle (a forward thrust propeller and coordination of the stern and 
rudder planes as ailerons). 

Here, attention is restricted to the longitudinal plane, and the formulation of section 4.1 


yield the state equations 


(m— X.)u  (Z, — m)wq +(Z, + mxs)q (UN = 12 E 
(m — Y, )vr 2 mzc(pr t d)  X,uuju + 
(B—W)s0 * X, (t) - X, 


Sy 


28 2 
X 5 190, + X,,5 UWwO, + Xu 0 


(ለያ) ያመ(ሙያ ሙለ ተእ/ጋተ የዘረ(1የመዘሆ-መ ሠሪ)” 
(mx; + M¿Jw + (2, — mx)uq +(mx — Y, )vp + 
4.23 
(A m~ L uw tb = eW aB xW cco ( 


M,w * M,q * MO) * M; wd. 
(m-Z,)w= (ms +Z,)g+Z,g+ Z,w+ (4, -m)vp+ 


(OE (IES Ip mz; p + q” )+ (4.24) 
(W - B)c8c6b + Z, (rt) + ጄ #6, 


88 





where X (t). Z,(1t) and M (t) are wave disturbance terms, and X,,, and Ö, are the forward thrust 


ihr 
and stern plane deflection angle, respectively. By making the further assumptions that 


d(1)=0 M,=M M, = M t Zy = Zu! Z7 ZUM 


ut ( HW Ww HW 


and eliminating all second order non-longitudinal plane terms yields the decoupled equations 


which will be used in subsequent discussions: 


(in — X;)u 2 (Z, — m)w0 + (Z, me )ፀ' ፦ mz 0 + 
X lulu + (B—W)s8+ X(t) + Xn, + (4.25) 


2ኗ 2 
X s uqO, t X,g wÒ, + X y 6, 


(1, — M¿JO = —m2z¿U— mz w8 + (mx; + M,)w+ 
(M,,, — Z, — xg )u8+(M,,, +X, — Z, uw + (4.26) 


(<,B—2,W)s0 + (x,B-xgW)cO+ M(t) + My WS, 


um 


(LW MA + Z; )ፀ +Z uw ዝሪ. MER: )u8 + 
T : (4.27) 
mo O OL O ZO. 


7=wcos8 — usin 0 (4.28) 
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Chapter 5 CONTROLLER DESIGN AND 
SIMULATIONS 


Shding control has been used successfully for position. and trajectory control of 
underwater vehicles by Yoerger and Slotine (1985, 1986) with the JASON remotely operated 
underwater vehicle, Healy and Marco (1992) with the Naval Postgraduate School AUVII, and 
Hills and Yoerger (1994) with the NUWC LDUUV. 

In this chapter, computer simulations are used to investigate the performance of 
variations of sliding control as applied to the pitch axis of the 21UUV. Robust sliding control 
routines are first developed for integrated depth, pitch, and forward speed control of the 21UUV, 
and then an adaptive sliding controller is introduced for the pitch axis. Coordinated control of the 
21UUV in the longitudinal plane is demonstrated by simulating the 21UUV making a depth 
change maneuver. The performance of the adaptive pitch controller in the presence of varying 
degrees of parametric uncertainty 1s also demonstrated. 

Disturbance cancellation properties of extensions to the adaptive pitch controller is then 
investigated, using first a monochromatic pitch disturbance of known frequency, and then a 
stochastic disturbance of known spectrum. 

The development of the robust sliding controller and the adaptive sliding controller 


presented here is found in (Slotine and Li 1991). 
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5.1 Robust Sliding Control 


5.1.1 Overview 


The single input dynamic system of the form 


x) = f(x)+b(x)u (5.1) 


is considered here, where the scalar x is the output, n is the order of the system and denotes the 
number of derivatives of x with respect to time, u is the scalar input, and the functions f and b are 


Siar 
x) and any other 


generally nonlinear functions of the state vector x=[x x 
measurable quantity. While the state vector x is assumed to be known exactly, the function f is 


not. Rather, the difference between f and its estimate is assumed to have a known bound that is a 


continuous function of the state and any other measurable quantity, 1.6. 


SX) (3.2) 





f- 


Similarly, it is assumed that the control gain b is not exactly known, but is of known sign and is 


bounded by known, continuous functions of x. 
The desired, realizable state trajectory is denoted x,, and the trajectory error vector is 


le 


denoted kx =x-x, =[x x . The hyperplane S(t) in the state space R'" is defined 


by 
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Te HEIC RP ET IU (5.3) 


where A 1s a strictly positive constant, and can be interpreted as the control bandwidth of the 
controller. 


Since b has known sign and bounds, and assuming here that b is positive, it can be seen 


that 
Üben De run (5D 
and that the bounds on / can be written 
b 
pos sp (6:5) 
where 
b=(b,,b.,. (5.6) 
and 
Bee, (5.7) 


While a control law which maintains s=0 would be ideal, 1t 1s also discontinuous across 
the hyperplane S(t) because of the uncertainties in the dynamic system. Thus, in practice, 
implementing such a law produces chattering in the control activity which is normally 
undesirable. To smooth such a discontinuous control law, a boundary layer neighboring the 


hyperplane can be used, namely 


B(t)={x, |s(x:t)|< @} 
92 0 


(15) 
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where the boundary layer thickness, Y can be made time varying to exploit the bandwidth of the 


system. 


It is shown in (Slotine and Li 1991) that the control law 


T bá RE sat(s/)| 


where 
ዘ = ay መ $4 xU? 
k =k(x)—k(x,) + A®/B(x,) 


and 


k(x) - p| F(x) +n]+(ß- Dál 


(5:9) 


(5710) 


(yl) 


(512 


with N a small, strictly positive constant, will provide robust stability in the presence of the 


assumed parametric uncertainty in f and b. The sat (saturation) function of equation (5.9) is 


given by 


Lx] =] = Sal ex 


a >] = sat(x) = sign(x) 


Boundary layer dynamics are given by 


Kk(x,)210/B(x,) = = P=-AD+ B(x, )k(x,) 
k(x,)<A@/B(x,) => b=-Ab/B?(x,)+k(x,)/B(x,) 


with 


@(0) = B(x, k(x, (0))/A 
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(9.1 >) 


(5.14) 


(5.15) 





5.1.2 Application to 21UUV Longitudinal Plane Equations 


In the derivation of control laws for the 21UUV, the state vectors Nand v of equation 
(4.1) are assumed to be measured and known exactly, while v is assumed unknown but of known 
bounds. 


Considering first 21} UUV depth control, from equation (4.28) 


Z=wcosd-usind -(wsind+ucosd) 8 (5216) 
es 
T b $, 
Then, 
፣ ጋጌ (3:17) 
P 
j (5.18) 
F, = |w| „cos 8 + |e] sin|| 
and 8. =1. 


From equation (5.16) it 1s seen that the depth control of this AUV ıs reliant upon its pitch 
control, which is an intuitive result given the actuators available on the vehicle. Therefore, the 
resulting control requirement for 0 affects not only the pitch angle, but the depth control of the 
vehicle as well. While 6 , 18 derived by the depth controller, ©, and O, are both quantities used 
in the derivation of the pitch controller. O, can be found by numerical integration of 0,. and O, 
can either be approximated by using a numerical derivative of 6 ¿» 9r can be uniformly replaced 
by 0 (with an accompanying decay in performance). Since 0 , 1$ à formulated quantity and is 
constructed using a smooth control law, taking the numerical derivative of 0 q 1$ a viable option 
for the derivation of the pitch controller. 


The control of the pitch angle 1s considered next, and from equation (4.26), 
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O= fy by, (5119) 


where 
fo= U,—M, m [ —mz 4 — mz," + (1x, + M,.)w+ 
ds zs mx, )u0 CM FA, ZZ iY (2220) 
(z,B — z;,W)sin0 - (x4B — x;W)cos6] 
b, -(,- Mj) Mg (5.21) 
and 
l = Ô, (322) 


The external disturbance due to waves, M_,(t), 1s assumed to be negligible for the moment. 
Then, where the estimates to the hydrodynamic and 21UUV body coefficients are 


denoted by a ^, 
= DUO (5.23) 


n E d. - M, Y [- m2 ¿we + (M,, = 2 nos ዘፀ Er 
(M,,, a5 X, E Z; uw ከ 2ቨ/) sin O + (324) 
(£,B - £,W)cos8] 


If the above hydrodynamic and body coefficients are known to within a certain range of values, 


their estimates can be taken as the algebraic mean of the highest and lowest values, and with 


|"... and EA known, F, can be found from equation (5.2). If this method ıs used to calculate 


ል 


by, equations (5.6) and (5.7) may not hold true, but by using equations (5.5) and (5.21), a 


conservative value for B, can be calculated, and is found to be a constant. 
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The control of the 21UUV forward speed is considered last, and from equation (4.25) 


ዘ= /, ተ (53.25) 
where 
ር በ ጋ ን (በሪ m)w8 + (Z; * mXg e - #ርራፀ + 
T (5.26) 
X „luju - (B — W)sin0 X ga 90, + X,5 UwO, + XS 02 
Don oe (5.27) 
and 
u,— X,, (528) 
Again, the external wave disturbance, X,,(1), is assumed to be negligible for the moment. 
Noting that equation (5.25) 1s of first order, 
sau (5.29) 
f,=(m-X,)'[(Z,, — m)w0 + Z + it, )0° + X, {eee + 
| (5.30) 


(B-W)sin0 +X,5 90, + X 5 UwO, + X 78, | 


and with 9 assumed known, F, can be found. Again, a conservative value of B, can be 
max 


calculated and 1s found to be a constant. 


With the completion of the above formulation, the choice of A¿ and A, remains. 


Guidance for the selection of sliding controller bandwidth 15 given in (Slotine and Li 1991) with 


concern for structural resonant modes, neglected time delays, and sampling rates addressed. 
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While unmodeled sternplane actuator dynamics may determine a reasonable choice for A,, the 
value of A. becomes dependent upon that of A¿. Because the outer loop depth controller 
develops the desired inner loop pitch rate and angle, A. <A,/4 1s used so that A, may be 


neglected when considering the equation of motion for z. 


5.2 Adaptive Sliding Control 


5.2.1 Application to 21UUV Pitch Equation 

Another approach to developing control laws in the presence of parametric uncertainty 
associated with the hydrodynamic and body coefficients of the equations in section 4.2 is to use a 
model based adaptive sliding controller. Where it was previously assumed in section 5.1.2 that 
the parameters in equations (4.25) through (4.27) are possibly time varying but bounded by 
upper and lower limits, the use of an adaptive sliding controller assumes no known bounds on the 
21UUV equation coefficients, but rather that the coefficients remain constant. 

Because pitch control is critical not only to the attitude but also to the depth control of 
the 21UUV, the application of adaptive sliding control will be presented here for the pitch 
equation. 


Equation (4.26) can be written 


a0 +a,w0 + aud +a,wu +assin9+a¿cos8+d,(1)+d,M,(1) =-408, (9.31) 


where 


on 





ug 


mu "nz, Z;, tmxg; M 





C ae a ae (b — a, = 
-ፅ., -M; - Ms 
ሪ..” ጎዲ, ፦ 47... =»W- 2.5 x-W — xB 
a y ዘ un ar G B A, = Ne Xp (5:52) 
- M; A - Ms 
| SEM. IZG. 
dy = = ee, 


Estimates of the values of the above hydrodynamic and body coefficients can be obtained 


through experimentation or the use of published experimental data such as that in (Hoerner 1965, 


1975) for similar shaped bodies. Then if [v and |... are known, D such that 1... e pa 


MAX 
also be conservatively estimated, where for the moment, M ,(f£) 1s considered negligible. 


Using the adaptive sliding pitch controller developed in the appendix yields the control 


and adaptation laws 


ô = uu “(kos = Ya) 


\ (5.33) 
a=-IY sı 
where 
y E ¿AO wÔ 10 wu  sinO cos J 
^ ^ ^ ^ ^ ^ ^ T 
(a Ma a iaa E] (5.34) 


pom 0 
SS Di Sas Do) 


and ką is such that 
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kD, > D (5.35) 


As noted above, the adaptation gain matrix, DL, is symmetric and. positive. definite. In 
practice, a diagonal matrix with strictly positive entries can be selected which allows for easy 
tuning of the adaptation gain for each parameter estimate, à,, individually. 

The adaptive sliding controller presented here, and referred to in later discussion, uses a 
time invariant boundary layer, ®. Slotine and Coetsee (1986) present similar results for the 
adaptive sliding controller with a time varying boundary layer, P(t). 

By comparing the laws for robust sliding pitch control and adaptive sliding pitch control, 
similarities can be seen. For the case when s remains inside its boundary layer, reformulating the 


robust sliding control law results in 


0, —u ( dig s- Yà) se: (5.36) 


un ^ 


Koen (1) 


which can be compared to equation (5.33). The on-line adjustment of Ag, using the laws given 


in section 5.1.1 for the robust shding controller allows for better exploitation of the control 


"bandwidth" available (Slotine and Li 1991). 


5.2.2 Adaptive Wave Disturbance Cancellation 
Notwithstanding the stochastic nature of ocean water waves and the development of 
chapters 2 and 3, first considered is the case where the 21UUV is operating under a sea 


dominated by regular waves of a single known frequency, œ, with unknown phase and amplitude. 


Then it 1s seen that the wave disturbance term in equation (5.31) can be written 


d M ,(t) = a, sin(or) + a, cos(@r) (5.37) 
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and by the addition of wave disturbance cancellation terms to Y and corresponding coefficient 


estimates to a in equations (5.33) and (5.34), namely 


Y-|6,- A6 wO  uO wu sind  cosO sino! coso | 
088) 


^ ^ ^ ^ ^ ^ ^ ^ ^ T. 
a= [a Qj, 04 a, a a Qa á; | 


the adaptive cancellation of the single frequency wave disturbance is possible. 

In an attempt to extend the above idea to better cancel the disturbance caused by random 
water waves, additional disturbance cancellation terms which span some portion of the spectrum 
of the stochastic disturbance can be added to Y. While robust stability and performance 
guarantees can no longer be provided as they were previously, a potential improvement in pitch 


performance seems likely under at least some conditions, and this idea will be investigated. 


5.3 Simulation Results 


Simulation results are presented here which demonstrate the performance of the two 
variations of sliding control presented in sections 3.1 and 5.2. In all cases, the simplified, 


coupled 21 UUV longitudinal plane equations ((4.25) through (4.28)) are used. 


5.3.1 Additional Modeling Considerations 

In the derivation of the longitudinal plane dynamic equations for the 21UUV, and 
subsequently during the controller design, no regard was given to actuator dynamics. For 
simulation purposes here, both the main thruster and sternplanes are modeled as first order 


systems which saturate, as depicted in figure 5.1. 
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Figure 5.1 Sternplane and Thruster Dynamics Model 


Another modeling consideration previously ignored is that the moment produced by the 
sternplanes grows with time until reaching its final value, M, uò.. The Wagner or "growth of 
lift’ function given in (Woods 1961) can be used to describe the lift dynamics. To develop 90% 


of the final lift, the sternplanes must travel approximately 6 chord lengths (Newman 1977). 


For the 0.1 meter chord length 21 UUV sternplane at 2 m/s, 90% lift is generated in about 
0.3 seconds. Equating this time to twice the time constant of a first order linear system results in 
an approximate bandwidth of 6.7 rad/sec. Since this bandwidth is somewhat greater than that 
assumed for the sternplane actuators, the lift dynamics are ignored for the purpose of these 
modeling considerations. 

Table 5.1 contains nominal values of the body, hydrodynamic, controller, and actuator 
coefficients used in the simulations. The values for the body and hydrodynamic coefficients were 
derived using strip theory and assuming the vehicle is neutrally buoyant. For hydrodynamic 


actuator constants (those with a ô, subscript) and values of the centers of buoyancy and gravity, 


the values provided by NUWC for the LDUUV are used. 
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Body Coefficients 


1619 kg 
15882 N 
15882 N 
7100 kg-m* 
2 Xx 10m 
0.005 m 


Om 


Om 


Hydrodynamic Coefficients 
849 ke 

-3066 kg-m 

-7255 kg-m* 

121.5 kg-m 

-355.1 kg 

121.5 kg-m 

-1619 kg 


-201.5 kg/m 


-770.4 kg 
-99.3 kg/m 
-41.5 kg 
-16.9 kg/m 
-55.3 kg 
-16.1 kg/m 


-36.7 kg/m 


Controller Constants Actuator Constants 


0. l 25 A eS 2 Ay tnt 0 


Ag a A, 0.5 0, max 30° ð, min 


Me , Tle . ኻን 0.01 2n max 650 N A min 





ID, diag[ 4000, 10, 1, 100, 1, 0.1] 


Table 5.1 Nominal Body, Hydrodynamic, Controller and Actuator Parameters 


Unless otherwise stated, it is assumed that the body coefficients of table 5.1 are known to 
within 5% error (estimated coefficients are within 5% of the nominal coefficients), and that the 
hydrodynamic coefficients are known to within 30% error. Actual values of the hydrodynamic 
and body coefficients used for simulation dynamics are those of the nominal values. Estimated 
values of the hydrodynamic coefficients used by the robust sliding controllers in their model of 
the 21UUV dynamics are 130% of the nominal values except for added mass terms 
(X, M,, M,,, Z,, and Z,) which are 70% of the nominal values. Similarly, estimated values of 
the body coefficients are 105% of the nominal values except for the inertia terms (m and /,) 


which are 95% of the nominal values. 
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Adaptive sliding controller initial estimates for à vary and are detailed with each 
associated simulation. 

Simulation results for the conditions investigated here have indicated the following 
estimates for bounds on maximum accelerations: 


pe[., «0.2 m/s, — [iu] 


max 


«0.2m/s, || «0.1deg/s' (5.39) 


Then from equations (5.32) and (5.39), it can be seen that Id, (t)... < D=0.08. 


5.3.2 Depth Trajectory Following 


Figures 5.2, 5.3, and 5.4 show simulation results of the 21UUV making a 10 meter depth 
change with the three integrated robust sliding controllers developed in section 3.1.2 controlling 
depth, pitch, and forward speed. Despite parametric uncertainty, depth trajectory following ıs 
good, and provides a baseline against which subsequent simulation results can be compared. 
Figure 5.5 depicts kgg (t) for the robust sliding pitch controller during the depth maneuver. 

By comparison, figures 5.6 and 5.7 show simulation results of the 21UUV making an 
equivalent depth maneuver with the same depth and speed controllers, but with the pitch control 
law from equation (5.33). With T=0, no parametric adaptation occurs, i.e. à—- à(0). The 
parameter estimates are calculated using the same body and hydrodynamic coefficient estimates 
used in the previous simulation. This pitch controller is tuned to have approximately the same 
feedback gain, k,, as that of the robust sliding controller demonstrated in figure 5.3. With 
ky — 17 ky, (t), depth and pitch control depicted in figures 5.6 and 5.7 closely match those of 
figures 5.2 and 5.3. Speed control, while also nearly identical to the previous case, is of 


secondary interest, and is not depicted in this nor subsequent simulations. 
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Figures 5.8 and 5.9 depict the same depth maneuver as in the previous cases with the 
pitch control law from equation (5.33), and with F =F, from table 5.1. All other conditions 
from the previous simulation remain the same. The resulting parameter adaptation is seen to 
have a positive effect upon the pitch and resulting depth control of the vehicle. Figure 5.10 
depicts the adaptation of the vector à shown with the parameter s, and its boundary layer ®,. 
Convergence of the parameter estimates to the values they estimate is not observed, nor is it 
necessarily expected. Narendra and Annaswamy (1989) and Slotine and Li (1991) both provide 
a thorough discussion concerning the condition of persistent excitation under which parameter 
estimates do converge to their target values. 

Next considered is the performance of the adaptive pitch controller in the presence of 
more substantial parametric uncertainty. It is now assumed that the hydrodynamic coefficients 
from table 5.1 are totally unknown, with the lone exception of M; , which is still assumed to be 
known to within 30% error. Dry body coefficients are assumed known to within 5% error as 
before. For the initial estimate of the adaptive controller parameter vector a, the previous 
estimated values of the body coefficients and M; are used, with the remainder of the 
hydrodynamic coefficient initially estimated to be 0. Actual values of the hydrodynamic and 
body coefficients used for simulation dynamics remain unchanged. 

Figures 5.11 and 5.12 show simulation results of the 21UUV attempting the 10 meter 
depth change with robust sliding depth and speed controllers, and the pitch control law from 
equation (5.33) with no adaptation. As illustrated, this control scheme is unstable. By 
comparison, figures 5.13 and 5.14 illustrate the beneficial effect of adaptation for the same initial 
conditions. Pitch control is shown to be stable, leading to good depth trajectory following. 
Figures 5.15 and 5.10 can be compared to see the effect of the added initial parametric 


uncertainty on the adaptation of the vector a. 
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Figure 5.2 21UUV Depth Trajectory with Robust Sliding Pitch Controller. Depicted (clockwise 
from top left) are (1) desired and simulated vehicle depth, (2) desired and simulated vertical 
velocity, (3) desired and simulated pitch rate, which is the control variable used in the depth 


equation, and (4) the generalized error parameter, s,, and its boundary layer, ®,. It is seen that 


despite substantial parametric uncertainty, good depth trajectory following is achieved. Note that 


s, remains well bounded by ®,, 


indicates that the resulting vehicle depth trajectory is maintained within the expected bounds of 
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which 1s a design feature of the robust sliding controller, and 


performance given the level of parametric uncertainty present. 
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Figure 5.3 21UUV Pitch Response During Depth Maneuver. 
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Here, the desired pitch rate 1s 


generated during the simulation by the depth controller, and the robust sliding pitch controller 


provides sternplane commands to follow the generated pitch trajectory. The sternplane is seen to 


remain within its saturation limits, and because of the high bandwidth of the sternplane actuator, 


sternplane angle follows the desired values closely. As in the case of the z controller, the error 


parameter, s, remains bounded by ®,. 
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Figure 5.4 21UUV Speed Response During Depth Maneuver. The forward speed controller here 
attempts to regulate the vehicle at 2 m/s during the change in vehicle depth. The steady state 
forward speed error seen at the end of the run is due to the parametric error in the estimated 


constants that the controller uses in its speed model. The effect of thruster dynamics are seen as 
delivered thrust varies from desired thrust. As in the previous two cases, s, remains bounded by 
9 


"E 
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k-theta-eff for robust sliding pitch controller 
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Figure 5.5 Robust Sliding Pitch Controller &,,,, (1). Though Ky, varies with time, it is seen that 


ko, * 1 forthe duration of this depth change maneuver. 
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Figure 5.6 21UUV Depth Trajectory with Adaptive Sliding Pitch Control Law, I = 0. Here, the 
pitch control law from equation (5.33) 1s used with no adaptation of parameter estimates. This 
pitch controller 1s very similar to the robust sliding controller used previously, differing only in 


that the feedback gain Kg, (1) from the robust sliding controller is approximated by a constant, 


k, = 7. Speed and depth controllers remain unchanged. As expected, depth trajectory following 


is very similar to that depicted in figure 5.2. 
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Figure 5.7 21UUV Pitch Response with Adaptive Sliding Pitch Control Law, I' 2 0. Pitch 


controller constants are k¿=7, ®,=0.012. a=a(0) is derived using the same parameter 
estimates used by the robust sliding controllers in the previous simulation. By design, pitch 


response 1s very similar to that depicted in figure 5.3. 
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Figure 5.8. 21UUV Depth Trajectory with Adaptive Sliding Pitch. Control Law, =I). Here, 


the pitch control law from equation (5.33) 1s used with adaptation of parameter estimates, and 


speed and depth controllers remain unchanged. Depth trajectory following is slightly improved 


from the previous two cases (figures 5.2 and 5.6) as a result of parameter adaptation by the 


adaptive pitch controller. 
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Figure 5.9 21UUV Pitch Response with Adaptive Sliding Pitch Control Law, F = F. Pitch 


controller values and à(0) remain the same as the previous simulation, but parameter adaptation 
now takes place. Pitch response is improved from the previous two cases (figures 5.3 and 5.7), 


resulting in the better depth trajectory following seen in the previous figure. 
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Figure 5.10 Adaptation of a and s,with ®,. The parameters å (t) through (t) are shown 


here for the previous simulation, and it is seen that adaptation occurs only while Isal > P}. It is 


noted that the parameter estimates, à;, do not necessarily converge to their target values, but s, 


does converge to the region inside boundary layer, P} which implies that 5,4— 0. 
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Figure 5.11 21UUV Depth Trajectory with Adaptive Sliding Pitch Control Law, I = 0. Here, 
the pitch control law from equation (5.33) 1s used with no adaptation of parameter estimates, and 
the depth and speed controllers remain unchanged. In this case, it is assumed that there is no 
prior knowledge of the hydrodynamic coefficients except for Ms , which is known to within 30% 
error as before. Body coefficients, as before, are assumed known to within 5% error. Because 


parametric uncertainty 1s so high, without adaptation this control scheme is unstable. 


114 





itch angle 









150 
100 


50 


theta and theta-d (--) (deg) 
q and q-d (--) (deg/s) 


50 100 


time (sec) 


stern plane angle 





200 
4 -200 E 
7 ፡ 
= -400 3 
I D 
= = 
= -600 J 
I 

-800 

0 50 100 
time (sec) 


80 pitch rate 


60 - ac fs 
40! ር 


20፦ ፡ = 








50 100 


time (sec) 


s and boundarv laver 





0 50 100 


time (sec) 


Figure 5.12 21UUV Pitch Response with Adaptive Sliding Pitch Control Law, I' 2 0. The pitch 


control law from equation (5.33) is used, and controller constants remain &, — 7 and €b, — 0.012. 


a = a(0) is derived using an assumption of no prior knowledge of hydrodynamic constants except 


for Ms , which is known to within 30% error as before. The estimates of body coefficients 


remain unchanged. Parametric uncertainty and the lack of parametric adaptation leads to 


instability of this pitch controller. 
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Figure 5.13. 21UUV Depth Trajectory with Adaptive Sliding Pitch Control Law, F = F. Here, 


the pitch control law from equation (5.33) is used with adaptation of parameter estimates. 


Controller constants remain unchanged, and à(0) is computed as in the previous simulation with 
no prior knowledge of hydrodynamic coefficients except for Mz , which is again known to within 


30% error. Despite high initial parametric uncertainty, this controller executes a depth change 


with good trajectory following. 
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Figure 5.14 21UUV Pitch Response with. Adaptive Sliding Pitch Control Law, I = E, Pitch 


controller values are the same as reported in figure 5.12, as is à(0). Despite high initial 
parametric uncertainty, this adaptive sliding pitch controller provides stable performance which 


translates to the good depth trajectory following seen in figure 5.13. 
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Figure 5.15 Adaptation of â and są with ®,. The parameters å (t) through 4¿(1) are shown 


here for the previous simulation with almost total hydrodynamic coefficient uncertainty. Like the 


case presented in figure 5.10, adaptation occurs only while 15,1 » ®,, and parameter estimates do 


not necessarily converge to their target values, but sy, > O. 
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5.3.3 Pitch Control in Regular Waves 

Now considered is the case where the objective of the pitch controller is to regulate the 
2] UUV pitch angle at O degrees in the presence of a single frequency disturbance, such as that 
provided by regular waves. Though the coupled equations of motion and forward speed control 
are aS before, the depth controller is disabled and no longer provides desired pitch information. 
Instead, 6,, 0,. and 0, are maintained at 0, while 4, is 2 m/s as before. This is done in an 
attempt to maintain the coupling effects of the speed and heave equations to the pitch equation, 
particularly in the presence of disturbances to all three axes, as will be investigated in section 
5.3.4. Such a situation of regulating pitch at O degrees may arise in practice during a phase of 
AUV operation when the stability of vehicle attitude takes priority over vehicle vertical position 


in the water column. 


In the following simulations, kọ and «b, are adjusted from the previous simulations to 
demonstrate the effect of the disturbance adaptation. The previous values of &, and «b, provided 
a pitch controller insensitive to the level of disturbance, which, when canceled, did not saturate 
the sternplane actuators at the 21 UUV forward speed of 2 nvs. 

Figure 5.16 depicts a simulation of the 21UUV subjected to a pitch disturbance of single 
frequency, with the pitch control law from equation (5.33) and no adaptation. In figure 5.17, 
pitch control of the 21UUV ıs shown with the same controller and subjected to the same 
disturbance, but with adaptation of the vehicle parameter estimates enabled. Though vehicle 
parameter estimates are updated through adaptation, little improvement in pitch regulation 1s seen 
with this second simulation. 

Assuming that the pitch disturbance frequency 1s known allows for use of the adaptive 
controller with disturbance cancellation terms as described in equation (5.38). Figure 5.18 
depicts this adaptive pitch controller performance in the presence of the same monochromatic 


disturbance used in the previous two simulations. The advantage of being able to include 
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cancellation terms which allow for phase and magnitude matching of the disturbance is seen in 
figure 5.18 as this controller substantially eliminates the effects of the regular waves and best 


maintains the 21 UUV pitch at 0 degrees. 
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Figure 5.16 Adaptive Pitch Control Law with Single Frequency Disturbance, F = 0. For this 
and the following simulations where a single frequency pitch disturbance is used, 4, = 80, 
o, — 0.001, and M, (f) 2 400 sin(1.7 t --0.3) N- m. The vector à(0) is calculated from here on 
with the assumption that body and hydrodynamic coefficients are known to within 5% and 30% 


error, respectively, with the estimated coefficient values chosen as described in section 5.3.1. 
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Figure 5.17 Adaptive Sliding Pitch Controller with Single Frequency Disturbance. This 
adaptive controller includes no disturbance cancellation terms, and [ = 1 ነ. While parameter 
estimate adaptation occurs, it is seen that pitch performance is very similar to the previous 


simulation with no adaptation. 
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Figure 5.18 Adaptation to a Single Frequency Disturbance by the Pitch Controller. Here, T = 
diag[T,, 1001]. Performance is improved with this controller which uses two disturbance 
cancellation terms, à, sin(1.7 ፤) and à, cos(1.7 r), to adapt to and cancel the pitch disturbance of 
known frequency, but of unknown phase and magnitude. Pitch response is significantly reduced 
from the previous cases, with final control action magnitude comparable to that of the previous 


two controllers. The error parameter, s, is seen to decrease in amplitude so that it becomes 


bounded by ®,, which implies that sy, — 0. 
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5.3.4 Pitch Control in Random Waves 

Considered here is pitch control for the 21UUV at a depth of 10 meters, in water 30 
meters deep, and in the presence of a disturbance which approximates that of sea state 2 
conditions. The theory presented in chapter 2 1s used to calculate the spectra of pitch, heave, and 
surge disturbances. A set of possible time realizations for the three disturbances is then 
generated. These disturbances are used in simulations to investigate the performance of the 
adaptive sliding pitch controllers presented in section 5.2. 

Figure 5.19 depicts the single-sided pitch, heave, and surge spectra. The first 60 seconds 
of pitch, heave, and surge time records used in the subsequent simulations are depicted in figure 
5.20. Associated with figure 5.20 are maximum horizontal and vertical displacement parameters 
of 0.2, justifying the use of the inertia dominated wave force assumption. 

Control simulations follow, and figure 5.21 shows an adaptive sliding pitch controller 
with no disturbance cancellation terms regulating 21 UUV pitch angle at O degrees in the presence 
of the stochastic wave disturbance. Attempting to partially reproduce the adaptive disturbance 
cancellation results seen in section 5.3.3, seven disturbance cancellation frequencies which span 
the known disturbance spectrum are added in feedforward, and their amplitudes and phase are 
adapted on-line. Figure 5.22 depicts this adaptive controller performance under the same 
conditions as those used with the previous simulation. The improvement in pitch regulation by 
this adaptive controller with disturbance cancellation terms is seen. Improvements for adaptive 


controllers with fewer disturbance cancellation terms, or for lower gain controllers with larger 


values of ®,, were found to be far less significant. 
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Figure 5.19 Wave Disturbance Spectra. The single-sided spectra above were generated using 
the ITTC recommended surface wave spectrum and the transfer functions for the 21UUV in 
inertia dominated wave forces presented in chapter 2. The dashed lines represent the head sea 
disturbance spectra for the AUV at O m/s, while the solid curves account for the shift of the 
spectra due to vehicle forward motion and frequency of encounter with the waves. As a spectrum 
is shifted to the right and its energy is spread over a larger range of frequencies, its peak 


amplitude drops, but the areas under the unshifted and shifted spectra remain equal. 
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Figure 5.20 Possible Time Realizations for Pitch, Heave, and Surge Disturbance. These 
disturbances were generated using the spectra in figure 5.19 and by superimposing 350 sine 
components as described in section 2.1.2. While the phase of each sine wave is random, proper 


relative phase between the sine components of each of the three disturbances (as presented in 


section 3.3.3.1) 1s preserved. 
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Figure 5.21 Adaptive Sliding Pitch Control with Stochastic Disturbance. Here, D = Dy 
k, = 200 and ®, = 0.001. Pitch angle regulation is good with this high gain controller, with l0l € 
0.2° during the simulation. This controller performance is used to measure against the following 
pitch controller which adds adaptive sine and cosine disturbance cancellation terms in an attempt 


to partially eliminate the effect of the random wave disturbance. 


126 





disturbance canceling a-hat 


| 0.5 
1111 lA 





theta and theta-d (--) (deg) 


-0.5 
-] 
0 50 100 0 50 100 
time (sec) time (sec) 
40 stern plane angle 4 x103 sand boundary layer 
'= | 
20 ርን | | | 
Bo Es ! £ à MATIN x |! 1111] = | | 
EMI T s "NE VL 
3 = | j- | | | j WI I | | 
-20 = -2H || 
-40 zz 
0 50 100 0 50 100 
ume (sec) time (sec) 


Figure 5.22 Adaptive Sliding Pitch Control with Stochastic Disturbance. As in figure 5.21, 
Kg = 200 and D, =0.001, but here 14 disturbance cancellation terms are added in feedforward 
and adapted on-line in an attempt to partially cancel the stochastic wave disturbance, and F = 
diag[1”, 1001]. The cancellation terms include a sine and cosine each at frequencies of 1.4, 1.5, 
1.6, 1.7, 1.8, 1.9, and 2.0 rad/sec. The coefficients of the 14 disturbance cancellation terms are 
plotted in the upper right. It is seen that steady state is not reached for the 14 coefficients, nor is 
it expected to be reached because exact cancellation of the stochastic disturbance with seven 
sinusoids is not possible. However, pitch amplitude is reduced over the previous simulation as 
0| € 0.17. The improvement in pitch regulation for adaptive controllers with fewer disturbance 


cancellation terms, or for lower gain controllers with a larger value of ®,, was not found to be 


nearly as pronounced. 
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The simulations demonstrating improvement in pitch regulation by adaptive sliding 
controllers with disturbance cancellation terms seen in the previous monochromatic and random 
disturbance cases each required higher gain controllers than the pitch controller used for depth 
trajectory following. In general, higher gain 1s less desirable because 1t makes a controller more 
susceptible to state measurement noise, and is more likely to excite unmodeled or non-ideal 
actuator dynamics. Particularly for the 21UUV at low forward speeds, it is apparent that a 
higher gain pitch controller would sooner cause sternplane saturation aS wave severity increased. 
Simulating the 21UUV operating at 5 meters depth in sea state 2 illustrates this point. 

Figure 5.23 depicts the 21UUV adaptive pitch controller (Ag =7, «b, —2 0.012, and 
F =T) regulating AUV pitch in the presence of a wave disturbance which simulates that of the 
vehicle operating in sea state 2, head seas, 30 meter deep water and at a depth of 5 meters. All 


three coupled longitudinal plane equations are again simulated so that the coupling effects of 


surge and heave on pitch are present. Figure 5.24 shows the adaptive pitch controller (A, = 200, 
co, = 0.001 and IT’ = diaglI p, 1001]) with 14 disturbance cancellation terms simulated under the 


Same conditions. As predicted, the higher gain causes actuator saturation, and the resulting 


controller performance is poorer than that of the lower gain pitch controller shown in figure 5.23. 
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Figure 5.23 Adaptive Sliding Pitch Control with Stochastic Disturbance. No feedforward 


cancellation terms are used in this controller, and I' 2 I, &, 27 and ®, =0.012. The simulated 
disturbance here is on the 5 meter deep 21UUV in sea state 2. Sternplane action is significantly 


less than in the following simulation which uses a higher gain pitch controller. 
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Figure 5.24. Adaptive Sliding Pitch Control with Stochastic Disturbance. This controller uses 14 
terms with 7 frequencies (1.4, 1.5, 1.6, 1.7, 1.8, 1.9, and 2.0 rad/sec) in an attempt to adapt to 
the stochastic disturbance. Here, I = diag[T',, 100I], X, 2200 and €, — 0.001 
saturation is evident and the resulting pitch performance is much poorer than that depicted in 


figure 5.23 for the lower gain controller. 
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5.3.5 Summary 

The sliding controllers developed in sections 5.1 and 5.2 for 21UUV longitudinal plane 
control were demonstrated here in simulation. Robust sliding control and adaptive sliding control 
techniques were applied to the pitch axis of the AUV, while robust sliding control laws alone 
were used for the depth and forward speed controllers. Development of the robust sliding control 
laws assumed that the hydrodynamic and body coefficients of the dynamics equations were 
bounded by some known, possibly time varying function, whereas the assumption implicit with 
the adaptive sliding control law was that the hydrodynamic and body coefficients remamed 
constant. 

21UUV pitch controller performance was shown to directly affect depth trajectory 
following as variations of sliding pitch control were used in the 3 axis, integrated control 
simulation of the 21UUV making a depth change maneuver. The effect of parametric uncertainty 
upon the performance of the pitch controllers was investigated, and the adaptive pitch controller 
was shown to provide good performance, resulting in good depth trajectory following, despite 
almost total hydrodynamic coefficient uncertainty. 

An extension of the adaptive pitch controller was shown to adapt to and cancel a pitch 
disturbance simulating that of a monochromatic wave of known frequency, but of unknown 
magnitude and phase. 

The developments of chapters 2 and 3 were used to generate a coordinated set of 
stochastic pitch, heave, and surge disturbances which the 21UUV might encounter in practice. 
The adaptive controller was again extended to include additional disturbance canceling terms in 
an attempt to achieve better pitch control performance. For the case presented, pitch regulation 
was improved by the addition of disturbance cancellation terms, though a high gain controller 


was needed to demonstrate this effect. 
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A second, more severe set of stochastic pitch, heave, and surge disturbances were then 
generated to demonstrate the effect of sternplane actuator saturation on the higher gaın adaptıve 
controller. When compared to the lower "bandwidth" pitch controller, the high gain controller 
was sooner to cause sternplane saturation with increasing disturbance severity, which resulted in 


its poorer performance when compared to the lower gain controller. 
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Chapter 6 CONCLUSIONS 


6.1 Summary 


Chapter 1 began with the motivation for this research work, namely to investigate the 
effect of water waves on a slender body, autonomous underwater vehicle operating in a region of 
the water column where inertia dominated wave forces dominate. The purpose of this research ıs 
not purely academic; Operation in coastal waters and sea state 2 by such a vehicle, the NUWC 
21UUV, is contemplated. The research objectives of this work were outlined and are consistent 
with a preliminary study of the problem: to develop a simple model for predicting wave 
disturbances; to perform tests which investigate the validity of the assumptions made in developing 
the model; and to use the model to demonstrate the effect of wave forces on the vehicle when 
controlled by a variation of 1ts proposed controller. 

Chapter 2 used existing linear wave theory, slender body strip theory, and linear time 
invariant systems theory to develop a method to predict the effect of wave disturbances on a 
stationary, slender body underwater vehicle. First, the prediction of monochromatic wave effects 
was addressed. These results were extended using a statistical description of waves and an 
assumption about the superposition of wave forces when operating in an inertia dominated wave 
force regime. Formulae for calculating the magnitude and phase of the transfer function from 
surface water waves to the disturbances on five of the six 21UUV axes were presented. Only the 


disturbance transfer function for the 21UUV roll axis did not lend itself to this method of 
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calculation, while a procedure for predicting surge, heave, sway, pitch and yaw axes disturbance 
transfer functions was shown. 

Having a statistical description of the sea surface, and using LTI system theory with the 
previously developed transfer functions allowed for development of a spectral description of the 
wave disturbances on the AUV. By using an earlier described method to generate a time 
simulation of the sea surface associated with a given sea spectrum, the generation of a simulated 
time record of wave disturbances on the AUV was shown. 

Chapter 3 described all aspects of the experimental tests which were conducted to 
investigate the validity of the transfer function model for wave forces on a stationary slender body 
vehicle. The experimental setup, consisting of a tow tank with motorized sliding carriage, the 
AUV tow body model, force and wave sensors, and data collection equipment were detailed. For a 
1:4.188 scale model of the 21UUV, it was shown that experiments conducted in the wave tank 
could preserve the wavelength to vehicle length ratios expected of the full scale vehicle in sea 
waves. For tests entailing the towing of the 21UUV model under waves, Froude number scaling 
was used and full scale equivalent forward speeds up to 1.5 m/s were investigated. The trials 
conducted during the course of testing were outlined, and they included investigating wave forces 
on the model 21UUV at a variety of model depths, model speeds, wave aspects, wave amplitudes, 
and over a range of wave frequencies. The processing of the raw data collected during the tests 
was also described. 

Presentation of experimental data along with the predicted results followed. Wave force 
and moment magnitude and phase on the AUV model were presented in a transfer function format 
for the tests conducted on the stationary model, and force and moment magnitude were presented 
for wave disturbances on the towed model. 

The stationary model tests largely validated the earlier developed method of predicting 


wave forces and moments on the body, both in magnitude and in phasé. A systematic discrepancy 
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in the prediction of the heave force on the model vehicle in head seas and at higher frequencies was 
noted. This discrepancy was seen in both heave force magnitude and phase. The transfer function 
description of forces and moments on the vehicle was supported by the data as 1t was shown that 
the hydrodynamic forces on the body were linearly related to the wave amplitude. 

The results of the data collected on the towed AUV model showed a deviation from the 
stationary body theory for the surge and heave forces. Towing the model AUV at the speeds 
investigated had minimal effect on the pitch transfer function magnitude. This result supported 
using the "stationary body" method of predicting the pitch disturbance on the forward moving AUV 
in subsequent simulations. 

Chapter 4 presented equations used to describe the motion of the 21UUV in 6 degrees of 
freedom. Coordinate systems, rigid body dynamics, hydrodynamic forces, and 21UUV body 
symmetry were addressed. The wave forces investigated earlier were included in these equations 
by superimposing them in an external disturbance vector. The 6 degree of freedom model was 
simplified to include only motion in the AUV longitudinal plane, resulting in coupled surge, heave, 
and pitch equations for the vehicle. 

Chapter 3 presented a method of applying sliding control techniques to the 21UUV 
longitudinal plane equations developed in the previous chapter. The effects of parametric 
uncertainty and pitch controller adaptation were demonstrated. It was found that even with almost 
no knowledge of AUV hydrodynamic coefficients, an adaptive sliding pitch controller was able to 
provide good pitch control performance throughout a commanded AUV depth change. The same 
control law without the benefit of parameter estimate adaptation proved unstable. This finding 
suggests that in practice, good 21UUV controller performance could still be achieved despite 
numerous in-field AUV configuration changes, and without the benefit of additional hydrodynamic 


coefficient analysis. 
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Pitch control in the presence of wave disturbances was also investigated, and 1t was shown 
that an adaptive sliding pitch controller could eliminate the effect of a monochromatic pitch 
disturbance of known frequency, but of unknown phase and magnitude. An extension of this idea 
to a stochastic pitch disturbance of known frequency spectra was investigated, and could only be 
shown to provide significant disturbance cancellation improvement for a higher gain controller than 
would normally be used. The negative effect of actuator saturation on pitch control performance 
was shown, supporting the use of a lower "bandwidth" pitch controller when controlling AUV pitch 


in waves. 


6.2 Future Directions 


The use of nonlinear sliding control with slender body underwater vehicles has been shown 
to be effective in practice (Healy and Marco 1992), and is currently being tested for a precursor to 
the 211UUV (Hills and Yoerger 1994). The work presented here concerning wave effects on an 
AUV is preliminary in nature, and the problem of AUV control in the presence of wave 
disturbances begs further study. While a method of calculating the disturbances caused by direct 
head seas and beam seas was presented, the operation of an AUV in waves of arbitrary aspect 1s of 
general interest. Similarly, computer simulations in this paper were limited to the longitudinal 
plane alone. Full 6 degree of freedom simulations of the 211UUV will be used to investigate the 
coupling effects of motion in the axes ignored in this paper. 

Though simulations can provide preliminary indications concerning the possible 
performance of an AUV's controller in the presence of waves, these simulations cannot realistically 
capture the full range of hydrodynamic effects and actuator and sensor performance found in 
practice. Ultimately, the richest method of studying the effect of wave disturbances on the 21UUV 


will be the full scale operation of the vehicle beneath the waves. 
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Appendix 
Below is contained the derivation of a robust adaptive sliding control law as applied to the 
21UUV pitch equation. This derivation is adapted from (Slotine and Li 1991). The resulting 


control and parameter adaptation laws are presented in section 5.2.1. 


a8 +a,w0+a,u0+a,wu +a,sin® +a, cos8 +d, (1) =-u ò, 


V=tas,,++a'T'a>0 a, >0 DEA 
1 uas a=a-a 
V-s,(-Ya-d(t)—-4/9,) - aT à d (D| <D 


Y =[0, wÔ uÔ wu sin 6056] 


a=[a a 4 a, a, a. | 


-u/ó,- Yà-ks, >  6,=ዘ”(ዴዴ -ኘ8) k 0 


V-s,(Yà-d,(t)-ks)*a T à ä=ä 


VAT" 4, Nasa) Spa =>  A=-TY se 


ee dia > 
Soakoso = Rose , + ንባ 
| 2 
Ya -So add (t) -KpBo|so | - ፍለ 


k®,2D > V<-ks, <0 
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lt is assumed that w and u are bounded. 


V20 V<0 = V is bounded => Sọ, a are bounded 


Sy bounded = 0,0 are bounded, and if € m d , are bounded (assumed). => 0.0 are bounded 


à bounded => a bounded 


for u » 0, 6, is bounded = Ə is bounded, and if 0, is bounded (assumed) — s, is bounded 


— V isbounded > V>0 by Barbalat's Lemma — 544 30 
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